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FOREWORD 


Tables of the U.S. Standard Atmosphere to heights in excess of 300 km have appeared 
in various editions over the past 18 years, the most recent being the U.S. Standard 
Atmosphere, 1975 (COESA, 1975) to heights of 1000 km. These publications have 
resulted from the work stimulated and guided by the Committee for Extension to 
the Standard Atmosphere (COESA). This committee has a membership of about 
30 organizations in government, industry, research institutions, and universities. The 
particular organizational representations in COESA resulted primarily from either 
the organization’s scientific or administrative needs, or the organization’s potential, 
for contributions to a standard atmosphere. In several cases, however, an organiza- 
tion’s association with COESA stemmed largely from the interests and talents of 
individual scientists employed by that organization. 

During the course of this work, COESA was led by three co-chairmen, each associated 
with one of the three principal sponsoring government agencies. These three agencies 
and the co-chairmen were: National Aeronautics and Space Administration (NASA), 
Maurice Dubin; National Oceanic and Atmospheric Administration (NOAA), Arnold 
R. Hull; and United States Air Force (USAF), K. S. W. Champion. 

The requirements and viewpoints of the many organizations participating in COESA 
have shaped the format of the U.S. Standard Atmosphere, 1975 (COESA, 1975). 
The scientific foundation for this Standard, however, is based upon the views of the 
international community of scientists including, to a large extent, the members of 
the Working Group of COESA. This Working Group, comprised of individual scien- 
tists representing nearly all of the 30 organizations in the parent committee, has 
based the U.S. Standard Atmosphere publication on current scientific and technical 
knowledge, rather than on policy concerns. 

The Working Group established five task groups for the preparation of the 1975 
Standard; some members contributed to more than one task group. Each of these 
task groups was charged with reviewing a particular category of available atmospheric 
data, and with preparing recommendations for possible revisions of the U.S. Stan- 
dard Atmosphere, 1962 (COESA, 1962). Three task groups (Task Groups I, II, and 
III) dealt with data involving the three overlapping height regions: 50 to 100 km, 80 
to 200 km, and 140 to 1000 km, respectively. Task Group IV was charged with 
combining the recommendations of Task Groups I, II, and III into an internally-con- 
sistent model. Task Group V reviewed the present-day knowledge of minor con- 
stituents. 



The membership of each of the five task groups is listed in the Foreword of the U.S. 
Standard Atmosphere, 1975 (COESA, 1975). This book also includes an extensive 
section of detailed tables of the properties of the Standard, which applies to a mean, 
midlatitude atmosphere for heights from 0 to 1000 km. The rationale for the model 
and an extensive discussion of the minor atmospheric constituents are also included. 

By contrast, Defining Constants „ Equations , and Abbreviated Tables of the 1975 
U.S. Standard Atmosphere emphasizes the definitions required to generate the tables 
of atmospheric properties of the Standard. This report also contains an abbreviated 
set of these tables, a set of graphs depicting the height profiles of these properties, 
and a discussion of two properties not considered in the 1975 Standard, that is, geo- 
potential pressure scale height and geopotential density scale height. 



CONTENTS 


Page 

ABSTRACT 1 

FOREWORD iii 

SYMBOLS vii 

INTRODUCTION 1 

INTERNATIONAL SYSTEM OF UNITS 2 

BASIC ASSUMPTIONS AND FORMULAS 2 

Adopted Constants 2 

Equilibrium Assumptions 10 

Gravity and Geopotential Altitude 12 

Mean Molecular Weight 16 

Molecular-scale Temperature versus Geopotential Altitude 

(0 to 84.3520 km') i8 

Kinetic Temperature versus Geometric Altitude (0 to 1000 km) 19 

COMPUTATIONAL EQUATIONS 22 

Pressure 23 

Number Density of Individual Species 24 

Mean Molecular Weight (Above 86 km) 29 

Pressure and Mass Density (Above 86 km) 30 

Total Number Density 30 

Mass Density 31 

Mole Volume 31 

Scale Heights 31 

Mean Air-particle Speed 35 

Mean Free Path 35 

Mean Collision Frequency 36 

Speed of Sound 37 

Dynamic Viscosity 38 


v 



CONTENTS (Continued) 


Page 

COMPUTATIONAL EQUATIONS (continued) 

Kinematic Viscosity 40 

Coefficient of Thermal Conductivity 40 

SUMMARY OF TABULAR VALUES OF ATMOSPHERIC PROPERTIES. ... 42 

Sea-level Values 42 

Conversion of Metric to English Units 43 

Tables of Atmospheric Properties 43 

ACKNOWLEDGMENTS 44 

REFERENCES 49 

APPENDIX A-BOUNDARY-VALUE NUMBER DENSITIES OF 

ATMOSPHERIC CONSTITUENTS 51 

APPENDIX B-A SEGMENT OF AN ELLIPSE TO EXPRESS TEMPERATURE 

VERSUS HEIGHT 55 

APPENDIX C— GEOPOTENTIAL AND MOLECULAR-SCALE 

TEMPERATURE 57 

APPENDIX D-THE CALCULATION OF A DYNAMIC MODEL FOR THE 

1975 U.S. STANDARD ATMOSPHERE 61 


vi 



SYMBOLS 


a One of the constant coefficients used to specify the elliptical segment of the 

temperature-height profile, T(Z) 

a. A set of species-dependent coefficients which, along with values of b., are used to 
define the set of height-dependent functions, D. 

A Another constant coefficient used to specify the elliptical segment of T(Z) 

A. A reaction-dependent rate coefficient used along with values B. and C^, to define 
k. 

j 

b A dimensionless subscript designating a set of integers { 0, 1 , 2, 3 . . .} with 0 

specifying sea-level conditions 

b. A set of species-dependent exponents which, along with values of a i? are used to 
define the set of height-dependent functions, D. 

B. A dimensionless, reaction-dependent exponent used in the expression for k. 

C. A reaction-dependent rate factor used as a part of an exponential expression in 
the definition of k. 

C $ The height-dependent speed of sound 

D. The set of height-dependent, species-dependent, molecular-diffusion coefficients 
for the atmospheric gas species O, 0 2 , Ar, He, and H 

f(Z) The hydrostatic term in the height-dependent expression for n. 

F. The dimensionless, sea-level, fractional-volume concentration of the ith member 

i 

of the set of atmospheric gas species 

F' Tire dimensionless, fractional-volume concentration of the i//t member of the set 

I 

of atmospheric gas species adjusted for 86-km height to account for the dissocia- 
tion of 0 2 


vii 



g The height-dependent acceleration of gravity (g(Z)) for 45° 

g^ The adopted constant involved in the definition of the standard geopotential 

meter, and in the relationship between geopotential height and geometric height 

H The geopotential height used as the argument for all tables up to 84.8520 km' 

(86 km) 

H p The height-dependent, local pressure scale height of the mixture of gases com- 

prising the atmosphere 

Hp The height-dependent, local geopotential pressure scale height of the mixture of 

gases comprising the atmosphere, and dependent upon the single variable, T M 

The height-dependent, local density scale height of the mixture of gases com- 
prising the atmosphere 

The height-dependent, local geopotential density scale height of the mixture of 
gases comprising the atmosphere 

i A dimensionless subscript designating the ith member of a set of gas species 

j A dimensionless subscript designating the j th member of a set of chemical reac- 

tions 

k The Boltzmann constant 

k. The reaction rate of the ]th chemical reaction 

k t The height-dependent coefficient of thermal conductivity 

K The height-dependent, eddy-diffusion (or turbulent-diffusion) coefficient 

L The height-dependent, mean free path 

L k b A set of gradients of T with respect to Z 

L m b A set of gradients of T M with respect to H 


M 


The height-dependent, mean molecular weight of the mixture of gases constituting 
the atmosphere 



M. The set of molecular weights of the several atmospheric gas species 

n. The set of height-dependent number densities of the several atmospheric gas 

species 

N Tire height-dependent, total number density of the mixture of neutral atmo- 

spheric gas particles 

N a The Avogadro constant 

P The height-dependent, total atmospheric pressure 

P. The partial pressure of the i th gas species 

q. One set of six adopted sets of species-dependent constants, that is, sets q., Q. , u., 

U., w., and W., all used in an empirical, species-dependent expression for the flux 
term v./(D. + K) 

Q. See q. 

r Q Tire adopted, effective earth’s radius used to compute g(Z) for 45° North latitude, 

and used for relating H and Z at that latitude 

R* The universal gas constant 

S The Sutherland constant, used in computing fi 

t The height-dependent Celsius temperature 

T The height-dependent, Kelvin kinetic temperature, defined as a function of Z 

for all heights above 86 km, and derived from T M for heights below 86 km 

T c The temperature coordinate of the center of the ellipse defining a portion of 

T(Z) 

T m The height-dependent, molecular-scale temperature defined as a function of H 

for all heights from sea-level to 86 km 

T^ The exospheric temperature 

u. See q. 

1 M 


IX 



U. See q i 

v. The flow velocity of the ith gas species 

v m The height-dependent mole volume 

V The height-dependent mean particle speed 

w. See q. 

W. See q. 

Z Geometric height used as the argument of all tables at heights above 86 km 

Z c The height coordinate of the center of the ellipse defining a portion of T(Z) 

a. The set of species-dependent, thermal-diffusion factors 

j 3 A constant used for computing p 

7 A constant representing the ratio of specific heat at constant pressure to the 

specific heat at constant volume used to define C s 

T The ratio g 0 /g Q 

e A dimensionless factor relating F. to Fj 

77 The height-dependent kinematic viscosity 

X A coefficient used to specify the exponential expression defining a portion of 

T(Z) 

p The height-dependent coefficient of dynamic viscosity 

v The height-dependent mean collision frequency 

£ A function of Z used in the exponential expression defining a portion of T(Z) 

p The height-dependent mass density of air 

o The effective mean collision diameter used in defining L and v 


x 



r A height-dependent coefficient representing the reduced height of the atomic 

hydrogen relative to a particular reference height, and used in the computation of 
n(H) (number density of hydrogen) 

0 The vertical flux of atomic hydrogen 

4> g The potential energy per unit mass of gravitational attraction 

<f> c The potential energy per unit mass associated with centrifugal force 


xi 
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INTRODUCTION 

The U.S. Standard Atmosphere , 1975 (COESA, 1975) is an idealized, steady-state represen- 
tation of the earth’s atmosphere from the surface of the earth to 1000-km altitude, as it is 
assumed to exist in a period of moderate solar activity. For heights from the surface to 
51 geopotential kilometers (km'), this Standard is identical to the U.S. Standard Atmo- 
sphere, 1962 (COESA, 1962), and is based on traditional definitions. These definitions, 
especially for heights below 20 km', were developed over a period of many years, and do 
not necessarily represent an average of the vast amount of atmospheric data available today 
from observations within that height region. For heights from 51 to 84.8520 km (that is, 
51.413 to 86 geometric kilometers), the 1975 Standard is consistent with averages of 
present-day atmospheric data, as represented by the traditional type of defining parame- 
ters. These include the linearly-segmented temperature-height profile, and the assumption 
of hydrostatic equilibrium, in which the air is treated as a homogeneous mixture of the 
several constituent gases. 

At greater heights, however, where dissociation and diffusion processes produce significant 
departures from homogeneity, the definitions governing the 1975 Standard are more 
sophisticated that those used at lower altitudes. In this high-altitude regime, the hydro- 
static equation, as applied to a mixed atmosphere, gives way to the more general equation 
for the vertical component of the flux for individual gas species (Colegrove et al., 1965; 
Keneshea and Zimmerman, 1970), which accounts for the relative change of composition 
with height. This flux equation simplifies to the hydrostatic equation for the special case 
when the atmospheric gases remain well mixed, as is the situation below 86 km. 



The temperature-height profile between 86 and 1000 km is not expressed as a series of 
linear functions, as at lower altitudes. Rather, it is defined in terms of four successive 
functions chosen not only to provide a reasonable approximation to observations, but 
also to yield a continuous first derivative with respect to altitude over the entire height 
regime. 

Observational data of various kinds provide the basis for independently determining vari- 
ous segments of this temperature-height profile. The observed temperatures at heights 
between 1 10 and 120 km were particularly important in imposing limits on the selection 
of the temperature-height function for that region. At the same time, the observed densi- 
ties at 1 50 km and above strongly influenced the selection of both the temperature and 
the vertical extent of the low-temperature isothermal layer immediately above 86 km. 

The mean temperatures derived from data sets associated with successive height regions 
were not necessarily continuous. In spite of this situation, it is necessary, for purposes of 
continuity and of mathematical reproducibility of the tables of this Standard, to express 
the temperature in a series of consecutive height functions from the surface to 1000 km. 
The expression for each successive function depends upon the end-point value of the pre- 
ceding function, as well as upon certain terms and coefficients peculiar to the related height 
interval. This total temperature-height profile applied to the fundamental continuity 
models (that is, the hydrostatic equation and the equation of motion), along with all the 
ancillary required constants, coefficients, and functions, defines the U.S. Standard Atmo- 
sphere, 1975 (COESA, 1975). The specification of this definition without any justification 
in terms of observed data is the purpose of this document. 

The definition of this Standard is completely consistent with the tables of two international 
standard atmospheres, that of the International Civil Aviation Organization (ICAO, 1964) 
defined up to 32-km altitude, and that of the International Organization for Standardiza- 
tion (ISO, 1973) defined up to 50 km. 

INTERNATIONAL SYSTEM OF UNITS 

The U.S. Standard Atmosphere , 1975 (COESA, 1975) is defined in terms of the Interna- 
tional System (S.I.) of Units (Mechtley, 1973). A list of the symbols, names, and the 
related quantities of the applicable basic and derived S.I. units, as well as of the three 
non-standard metric units and one English unit employed in this Standard Atmosphere, 
is presented in table 1 . 

BASIC ASSUMPTIONS AND FORMULAS 
Adopted Constants 

For purposes of computation, it is necessary to establish numerical values for various con- 
stants appropriate to the earth’s atmosphere. The adopted constants are grouped into 
three categories. Category I includes those constants which are common to many branches 
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Table 1 

List of Units Applicable to the U.S. Standard Atmosphere, 1975 

(COESA, 1975) 


Symbol 

Name 

Quantity 

Basic S.I. 



m 

meter 

Length 

kg 

kilogram 

Mass 

s 

second 

Time 

K 

Kelvin 

Thermodynamic temperature 

mol 

mole 

The amount of a substance 

Derived S.I. 



N 

newton 

Force (kg • m/s) 

Pa 

pascal 

Pressure (N/m 2 ) 

J 

joule 

Work, energy or quantity of 
heat (N • m) 

W 

watt 

Rate of energy (or heat) 
transfer (J/s) 

Non-standard 



mbar 

millibar 

Pressure 1 00 (N/m 2 ) 

torr at 0°C 

torr 

Pressure 133.322 (N/m 2 ) 

°C 

Celsius 

degree 

Temperature (Kelvin minus 
273.15) 

m' 

geopotential 

meter 

Variable unit of length 

English 



ft 

foot 

Length (0.3048 m)* 


* Exact definition 


of the physical and chemical sciences, and are here considered to be fundamental constants. 
Some of these may be multi-valued as in the case of M; representing the molecular weight 
of the i th gas species. Category I includes three single-valued and one multi-valued constant. 
Category II includes those constants which, in addition to the Category I constants and a 
suitable set of equations, are sufficient to define that portion of the 1975 Standard Atmo- 
sphere below 86 km. This category includes nine single-valued and three multi-valued 
constants. Category III includes all the remaining constants which, along with Category I 
and Category II constants and the related equations required for the low-altitude portion, 
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plus an expansion of that set of equations, are necessary to define that portion of the 1975 
Standard Atmosphere above 86 km. This category includes 7 single-valued and 12 multi- 
valued constants. 

The various single-valued and multi-valued constants are listed alphabetically by symbol with- 
in each of the three categories. The discussion of each of these constants includes the nu- 
merical value and dimensions, except that in the case of multi-valued constants, the values 
are given in one or another of five tables immediately following the listing of the three 
categories of constants. 

Category I Constants 

k The Boltzmann constant, k = 1 .380 622 X 10“ 23 N • m/K, is theoretically equal 

to the ratio R*/N A , and has a value consistent with the carbon-1 2 scale, as cited 
by Mechtly (1973). 

M. The quantity M. represents the set of the 10 values (kg/kmol) of molecular 

weight of particular atmospheric gases, listed in table 2. This set is based upon 
the carbon-12 isotope scale for which C 12 = 12. This scale was adopted in 1961 
at the Montreal meeting of the International Union of Pure and Applied Chemis- 
try. 

N a The Avogadro constant, N A = 6.022 169 X 10 26 kmol" 1 , is consistent with the 

carbon-12 scale, and is the value cited by Mechtly (1973). 

R* The gas constant, R* = 8.314 32 X 10 3 N • m/(kmol • K), is consistent with the 

carbon-12 scale, but is not the value cited by Mechtly (1973). Rather, it is the 
value used in the U.S. Standard Atmosphere , 1962 (COESA, 1962). This value 
is not exactly consistent with the above-listed values of k and N A . 

Category II Constants 

F The general member, F., of the set of values of fractional-volume concentrations 

of atmospheric gases (table 2) represents the assumed relative concentrations of 
the several gas species comprising dry air at sea level. The values of this set are 
identical to those given in the 1962 Standard. Except for minor modifications 
which are based upon C0 2 measurements (Keeling, 1960), these values are the 
same as those given by Glueckauf (1951) and are based upon the earlier work of 
Paneth (1939). 

g 0 The quantity, g 0 = 9.80665 m/s 2 , represents the adopted sea-level value of the 

acceleration of gravity for this Standard. This value is the one originally adopted 
by the International Committee on Weights and Measures in 1901 for 45° lati- 
tude. Even though it has since been shown to be too high by about five parts in 
1 0,000 (List, 1 968), this value has persisted in meteorology and in some stan- 
dard atmospheres as the value associated with 45° latitude, even though it 
applies more precisely to a latitude of 45°32'33". 
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g 0 The dimensional constant, = 9.80665 m 2 /(s 2 • m'), selected to relate the 

standard geopotential meter (in') to geometric height is numerically equal to 
g 0 , but with appropriately different dimensions, that is, m 2 /(s 2 • m'). This 
constant implicitly defines 1 m' as the vertical increment through which one 
must lift 1 kg to increase its potential energy by 9.80665 joules (J). The geome- 
tric length of this vertical increment varies inversely with the height-dependent 
value of g. 

H b The general member, H b , of the set of geopotential-height values listed in 

table 3 represents the base of each of eight successive atmospheric layers, 
defined to have particular constant values of L M b . 

L m b The general member, L M b , of the set of seven values of gradient dT M /dH of 
molecular scale temperature T M (Minzner and Ripley, 1956) with respect to 
geopotential H listed in table 3 represents the fixed value appropriate throughout 
its related layer, H b to H b+1 . The corresponding pairs of values of H b and L M b 
are based partly on tradition and partly on present-day observations. The first 
five of these pairs are identical to those of the first five layers of the 1962 Stan- 
dard, while the remaining two pairs have been newly-selected to provide a rea- 
sonable fit to the presently-available atmospheric data. The first two values of 
the related sets have their origin in one of the earliest aeronautical standard 
atmospheres (Toussaint, 1919) and were approximated in the first U.S. Stan- 
dard Atmosphere (Diehl, 1925). 

P Q The standard sea-level atmospheric pressure, P Q = 1.013250 X 10 5 Pa (or N/m 2 ), 

was adopted in 1947 in Resolution 164 of the International Meteorological 
Organization. This corresponds to the pressure exerted by a column of mercury 
0.760 m high, having a density of 1.35951 X 10 4 kg/m 3 , and subject to an accel- 
eration due to gravity of 9.80665 m/s 2 . This equivalency definition was adopted 
by the International Commission on Weights and Measures in 1948. 

r Q The effective earth’s radius for purposes of calculating geopotential altitude at 

any latitude is readily obtained from equations given by Harrison (1968). The 
value of r Q = 6.356766 X 10 6 m used in this Standard corresponds to the lati- 
tude for which g = 9.80665 m/s 2 . 

T 0 The standard sea-level temperature, T Q = 288.15 K, is based upon two inter- 

national agreements. The first of these is Resolution 192 of the International 
Commission for Air Navigation which in 1924 adopted 15°C as the sea-level 
temperature of the International Standard Atmosphere. This value has been 
retained unchanged in all known standard atmospheres since that date. The 
second agreement is that of the 1954 Tenth General Conference on Weights 
and Measures which set the fixed point of the Kelvin temperature scale at the 
triple-point temperature 273.16 K, which is 0.01 K above the ice-point tempera- 
ture at standard sea-level pressure. 
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S The Sutherland constant, S - 1 10 K, (Hilsenrath et ah, 1955) is a constant in 

the empirical expression for dynamic viscosity. 

j3 The quantity 0 = 1.458 X 10 6 kg/(s • m • K I/2 ) (Hilsenrath et ah, 1955) is a con- 

stant in the expression for dynamic viscosity. 

7 The ratio of specific heat of air at constant pressure to the specific heat of air 

at constant volume is a dimensionless quantity with an adopted value y = 1.4. 
This is the value adopted by the Aerological Commission of the International 
Meteorological Organization in Toronto, 1948. 

o The mean effective collision diameter, o = 3.65 X 10' 10 m, of gas molecules is a 

quantity which varies with gas species and temperature. The adopted value is 
assumed to apply in a dry, sea-level atmosphere. Above 85 km, the validity of 
the adopted value decreases with increasing altitude (Herschfelder et ah, 1964; 
Chapman and Cowling, 1960) due to the change in atmospheric composition. 

For this reason, the number of significant figures in tabulations of quantities 
involving o is reduced from that used for other tabulated quantities at heights 
above 86 km. 

Category III Constants 

a 4 The quantity a. represents the general member of a set of five values (m' ! • s' 1 ) 

of a species-dependent coefficient listed in table 4, and used in equation 8 for 
designating the height-dependent, molecular-diffusion coefficient D . for the 
related gas species. (Seeh.) 

b i The quantity h represents the general member of a set of five values (dimension- 

less) of a species-dependent exponent listed in table 4, and used, along with the 
corresponding value of a i? in equation 8 for designating the height-dependent, 
molecular-diffusion coefficient for the related gas species. The particular values 
of a 4 and b. adopted for this Standard have been selected to yield a height varia- 
tion of D. consistent with observed number densities. 

K ? The quantity K ? = 1.2 X 10 2 m 2 /s is the adopted value of the eddy-diffusion 

coefficient K, at Z ? = 86 km, and in the height interval from 86 up to 91 km. 
Beginning at 91 km and extending up to 115 km, the value of K is defined by 
equation 7b. At 1 1 5 km the value of K equals K 10 . 

K J0 The quantity K J0 = 0.0 m 2 /s is the adopted value of the eddy-diffusion coeffi- 
cient K at Z 1Q = 120 and throughout the height interval from 1 1 5 km through 
1000 km. 
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L^ b The two-valued set of gradients L Kb = dT/dZ listed in table 5 was specifically 
selected for this Standard to represent available observations. The first of these 
two values of L„ . is associated with the layer 86 to 91 km, and the second with 
the layer 1 1 0 to 120 km. 

n(0) ? The quantity n(0) 7 = 8.6 X 10 16 m 3 is the number density of atomic oxygen 
assumed for this Standard to exist at Z 7 = 86 km. This value of atomic oxygen 
number density, along with other defined constants, leads to particular values of 
number density for N 2 , 0 2 , Ar, and He at 86 km. (See Appendix A.) 

n(H) n The quantity n(H) n = 8.0 X 10 10 m' 3 is the assumed number density of atomic 
hydrogen at height Z n = 500 km, and is used as the reference value in computing 
the height profile of atomic hydrogen between 150 and 1000 km. 

q. The quantity q s represents the first set of six sets of species-dependent coeffi- 

cients listed in table 6 (that is, sets of q., Q j; a, U, w., and W.) the correspond- 
ing members of all six of which are simultaneously used in an empirical expres- 
sion for the vertical transport term Vj/fDj + K) in the vertical flux equation for 
the particular gas species. The species-dependent values of all six sets have been 
selected for this Standard to adjust number-density profiles of the related gas 
species to particular boundary conditions at 1 50 and 450 km, as well as at 97 km 
in the case of atomic oxygen. These boundary conditions all represent observed 
or assumed average conditions. (See equation 37.) 

Q. The quantity Q s represents the second set of the six sets of constants described 

along with q. above, and listed in table 6. 

T 9 The quantity T 9 = 240 K represents the kinetic temperature at Z 9 = 1 10 km. 

This temperature has been adopted along with the gradient L K 9 = 12 K/km to 
generate a linear segment of T(Z) for this Standard between 1 1 0 and 1 20 km. 

T^ The quantity T^ = 1000 K represents the exospheric temperature, that is, the 

asymptote which the exponential function, representing T(Z) above 1 20 km, 
closely approaches at heights above about 500 km, where the mean free path 
exceeds the scale height. The value of T^ adopted for this Standard is assumed 
to represent mean solar conditions. 

u. The quantity u. represents the third set of the six sets of constants described 

along with q. above, and listed in table 6. 

U. The quantity U represents the fourth set of the six sets of constants described 

along with q. above, and listed in table 6. 
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w. The quantity w. represents the fifth set of the six sets of constants described 

along with q. above, and listed in table 6. 

W. The quantity W. represents the sixth set of the six sets of constants described 

along with q ; above, and listed in table 6. 

Z b The quantity Z b represents a set of six values of Z for b equal to 7 through 12. 

The values Z 7 , Z 8 , Z g , and Z ]0 correspond respectively, to the base of succes- 
sive layers characterized by successive segments of the adopted temperature- 
height function for this Standard. The fifth value, Z n , is the reference height 
for the atomic hydrogen calculation, while the sixth value, Z ]2 , represents the 
top of the region for which the tabular values of the Standard are given. These 
six values of Z b , along with the designation of the type of temperature-height 
function associated with the first four of these values, plus the related value of 
L„ . , for the two segments having a linear temperature-height function, are 
listed in table 5. 

a. The quantity a. represents a set of six adopted species-dependent, thermal-diffu- 

sion coefficients listed in table 4. 

4> The quantity (p = 7.2 X 1 0 11 m' 2 • s" 1 for the vertical flux is chosen as a com- 

promise between the classical Jeans’ escape flux for T^ = 1000 K, with correc- 
tions to take into account deviations from a Maxwellian velocity distribution at 
the critical level (Brinkman, 1971) and the effects of charge exchange with H + 
and 0 + in the plasmasphere (Tinsley, 1973). 


Table 2 

Molecular Weights and Assumed Fractional-volume Composition of Sea-level Dry Air 


Gas Species 

Molecular Weight 

Fractional Volume 

(kg/kmol) 

F. (dimensionless) 

N 2 

28.013 4 

0.780 84 


31.998 8 

0.209 476 

At 

39.948 

0.009 34 

co 2 

44.009 95 

0.000 3 1 4 

Ne 

20.183 

0.000 018 18 

He 

4.002 6 

0.000 005 24 

Kr 

83.80 

0.000 001 14 

Xe 

131.30 

0.000 000 087 

ch 4 

16.043 03 

0.000 002 

h 2 

2.015 94 

0.000 000 5 
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Table 3 

The Adopted Reference Levels and Gradients of the Linearly-segmented 
Temperature-height Profile from the Surface to 86 Geometric 
Kilometers and the Related Derived Values of 
Molecular-scale Temperature 


Subscript 

Ci copotential 
Height 

Molecular-scale 

Temperature 

Gradient 

Form of Function 
Relating T to 11 

Derived Values of 
Molecular-scale 
Temperature 

b 

H h tkm') 

L M h <K/km’) 


T M .h< K > 

0 

0 

-6.5 

Linear 

288.150 

1 

1 1 

0.0 

Linear 

216.650 

i 

20 

+ 1.0 

Linear 

216.650 

3 

32 

+2.8 

Linear 

228.650 

4 

47 

0.0 

Linear 

270.650 

5 

51 

-2.8 

Linear 

270.650 

6 

71 

-2.0 

Linear 

214.650 

7 

84.8520 



186.946 


Note: These values of ll b and l. M .h plus Ihc defined sea-level value of T (equal to Tm, 0> completely specify the gcopotcntial-heigh. profile of T M f»t» 
the surface to 86 geometric kilometers, as designated by the derived values of * n ,h ' s ta ^ c 


Table 4 

A Set of Species-dependent, Thermal-diffusion Factors and Two Other Sets of 
Species-dependent Constants Required in Specifying the Height-dependent 
Function of the Molecular-diffusion Coefficient for the 
Several Species Listed 


C i as 

a (dimensionless) 

ii. = (m 1 * s 1 ) 

b. (dimensionless) 

N 

0.00 



2 

o 

0.00 

6.986 X 10 20 

0.750 

o 

0.00 

4.863 X 10 20 

0.750 

2 

Ar 

0.00 

4.487 X 10 20 

0.870 

He 

-0.40 

1.700 X 10 21 

0.691 

H 

-0.25 

3.305 X 10 21 

0.500 
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Table 5 

The Constants and Functions Adopted to Define the Four-layer Temperature-height 
Profile for Altitudes Between 86 and 1000 km, Plus the Derived Temperatures at the 
Boundary Heights of the Several Layers 


Subscript 

( iconic trie 
l {eight 

Kinetic- temperature 
Gradient 

Form of Function 
Relating T to X 

Derived Kinetic 
Temperature 

b 

Z b (km) 

L k b (K/km) 

T (K) 

7 

86 

0.0 

Li n e ci r 

186.87 

8 

91 


I lliptical 

186.87 

9 

1 10 

12.0 

Linear 


10 

120 


Lxponential 

360.00 

1 1 

500 



999.24 

12 

1000 



1 000.00 


Note: These adopted specifications, including the adopted kinetic temperature, 240 K at 110 km, and the kinetic temperature, 186.87 K derived fromT\j 
at 86 km, plus the specific form of the exponential function, equation 31 . and the requirement that dT/dZ be continuous from Z ~ 86 km to Z - 
1000 km. define the height profile of T between these height limits. The specific form of the ellipse, equation 27, which satisfies the several adopted 
conditions is derived in Appendix IL 


Table 6 

Values of Six Sets of Species-dependent Coefficients Applicable to t lie Empirical 
Expression Representing the Flux Term v./(D. + K) in the Equation for 
Number Density of the Four Species Listed 


Gas 

q. (km‘ 3 ) 

Q, (km' 3 ) 

u. (km) 

U. (km) 

w. (km' 3 ) 

W. (km' 3 ) 

O 

-3.416248 X 10' 3 * 

-5.809644 X 10"* 

97.0 

56.9031 1 

5.008765 X 10 4 

2.706240 X 10~ 5 


0 

1.366212 X 10^ 


86.000 


8.333333 X 10* 5 

Ar 

0 

9.434079 X 10' 5 


86.000 


8.333333 X \Qr 5 

He 

0 

-2.457369 X I0 4 


86.000 


6.666667 X 10"* 


♦ This value of q; applies only for 86 < L < 97 km. For Z > 97 km, qj - 0.0 kin'L 


Equilibrium Assumptions 

The air is assumed to be dry, and at heights sufficiently below 86 km, the atmosphere is 
assumed to be homogeneously mixed with a relative-volume composition leading to a con- 
stant mean molecular weight, M. In this height region of complete mixing, the air is 
treated as if it were a perfect gas, and the total pressure P, the temperature T, and the total 
density p at any point in the atmosphere are related by the equation of state, that is, the 
perfect gas law, one form of which is 


p • R* • T 

P 

M 


( 1 ) 
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where R* is the universal gas constant. An alternate form of the equation of state, (equa- 
tion 1) this one in terms of the total number density, N, and the Avogadro constant, N A , is 


P = 



T 


( 2 ) 


This form represents the summation of P ( , the partial pressures of the individual gas species, 
where P is related to n., the number density of the i th gas species in the following expres- 
sion: 

Pi = n. • k • T (3) 


where k is the Boltzmann constant, which is equivalent to R*/N A . 

Within the height region of complete mixing, the atmosphere is assumed to be in hydrosta- 
tic equilibrium, and to be horizontally stratified so that dP, the differential of pressure, is 
related to dZ, the differential of geometric height, by the relationship 

dP = -g • p • dZ (4) 

where g is the height-dependent acceleration of gravity. The elimination of p between 
equations 1 and 4 yields another well-known form of the hydrostatic equation, which 
serves as the basis for the low-altitude pressure calculation: 

dP -g • M 

dCnP = — = — • dZ (5) 

P R* • T 

Above 86 km the hydrostatic equilibrium of the atmosphere gradually breaks down as dif- 
fusion and vertical transport of individual gas species lead to the need for a dynamically- 
oriented model including diffusive separation. Under these conditions it is convenient to 
express the height variations in the atmospheric number density in terms of the vertical 
component of the flux of the molecules of individual gas species (Colegrove et al., 1965). 
In terms of the i th gas species, this expression is 


iij • Vj + Dj 


where 


dn. n. • ( 1 + a.) 

i i v i' 

+ 

dZ T 


dT_ g ' n i ' M i \ 

dZ R* • T / 


/ ^ n i 

n. 

i 

dT 

g • n. * M \ 

\ dZ 

i 

T 

dZ 

f R* * T / 


( 6 ) 


0 


v. = the vertical transport velocity of the i th species, 

D. = the height-dependent, molecular-diffusion coefficient of the i th species diffusing 
through N 2 , 

a. = the thermal-diffusion factor (dimensionless) of the i th species. 


11 



M ; = the molecular weight of the i th species, 

M = the molecular weight of the gas through which the i th species is diffusing, and 

K = the height-dependent, eddy-diffusion coefficient. 

The function K is defined differently in each of three height regions: 

1. For 86 < Z < 95 km, 


K = K ? = 1.2 X 10 2 m 2 /s 


2. For 95 < Z < 115 km, 

K = K ? • exp 

3. For 115 <Z< 1000, 


400 


400 - (Z - 95 y 


(7a) 


(7b) 


K = K 10 = 0.0 m 2 /s 

The function D ( is defined by 


D. 

1 


a. 


i 


2n ; 



(7c) 


(8) 


where a i and are the species-dependent constants defined in table 4, while T and Zn. are 
both altitude-dependent quantities which are specified in detail below. The values of D. , 
determined from these altitude-dependent quantities, and the defined constants a. and b. 
are plotted in figure 1 as a function of altitude for each of four species, O, 0 2 , Ar, and He. 



MOLECULAR-DIFFUSION AND EDDY-DIFFUSION COEFFICIENTS, m ? /s 


Figure 1 . Molecular-diffusion and eddy-diffusion coefficients versus geometric altitude. 
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The value of D i for atomic hydrogen, H, for heights just below 150 km is also shown in 
figure 1 . This same figure contains a graph of K as a function of altitude. It is apparent 
that for heights sufficiently below 90 km, values of D. are negligible compared with K, 
while above 1 1 5 km, the reverse is true. In addition, it is known that the flux velocity, v. , 
for the various species becomes negligibly small at altitudes sufficiently below 90 km. 

The information regarding the relative magnitudes of v., D., and K permits us to consider 
the application of equation 6 in each of several regimes. One of these regimes is for 
heights sufficiently below 90 km, such that v. and D. are both extremely small compared 
with K. Under these conditions, equation 6 reduces to the following form of the hydro- 
static equation: 


dn 


i [ dT 
n. T 

i 


-g • M 
R* • T 


• dZ 


(9) 


Since the left-hand side of this equation is seen through equation 3 to be equal to dR/P., 
equation 9 is seen to be the single-gas equivalent to equation 5. Consequently, while equa- 
tion 6 was designed to describe the assumed equilibrium conditions of individual gases 
above 86 km, it is apparent that it also describes such conditions below that altitude. Here 
the partial pressure of each gas comprising the total pressure varies in accordance with the 
mean molecular weight of the mixture, as well as in accordance with the temperature and 
the acceleration of gravity. Nevertheless, equation 5, expressing total pressure, represents 
a convenient step in the development of equations for computing total pressure versus geo- 
metric height, when suitable functions are introduced to account for the altitude variation 
in T, M, and g. 


It has been customary in standard-atmosphere calculations to effectively eliminate the 
variable portion of the acceleration of gravity from equation 5 by the transformation of 
the independent variable Z to geopotential altitude H. This simplifies both the integration 
of equation 5 and the resulting expression for computing pressure. The relationship be- 
tween geometric and geopotential altitude depends upon the concept of gravity. 


Gravity and Geopotential Altitude 

Viewed in the ordinary manner from a frame of reference fixed in the earth, the atmosphere 
is subject to the force of gravity. The force of gravity is the resultant (vector sum) of two 
forces: (1) the gravitational attraction in accordance with Newton’s universal law of gravi- 
tation, and (2) the centrifugal force, which results from the choice of an earthbound, rota- 
ting frame of reference. 

The gravity field, being a conservative field, can be derived conveniently from the gravity 
potential energy per unit mass, that is, from the geopotential <h. This is given by 

4> = 4>g + <b c 0°) 
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where <J> G is the potential energy per unit mass of gravitational attraction, and 4> c is the 
potential energy per unit mass associated with the centrifugal force. The gravity per unit 
mass is 


g = V4> O') 

where V4> is the gradient (ascendant) of the geopotential. The acceleration due to gravity 
is denoted by g and is defined as the magnitude of g, that is, 

g=.|gj. = |V*| (12) 

When moving along an external normal from any point on the surface < t> l to a point on the 
surface <i > 2 infinitely close to the first surface, so that <f> 2 = + ^ ie incremental work 

performed by shifting a unit mass from the first surface to the second will be 

d<& = g • dZ 03) 

hence 

/ 

<i> = C Z g -dZ. (14) 

The unit of measurement of geopotential (Appendix C) is the standard geopotential meter 
which represents the work done by lifting a unit mass one geometric meter, through a 
region in which the acceleration of gravity is uniformly 9.80665 m/s 2 . 

The geopotential of any point with respect to mean sea level (assumed zero potential), 
expressed in geopotential meters, is called geopotential altitude. Therefore, geopotential 
altitude, H, is given by 


H _ <t> _ 1 

So So 



g • dZ 


(15) 


and is expressed in geopotential meters (m') when the unit geopotential, g^ , is set equal to 
9.80665 m 2 /(s 2 • m'). 


With geopotential altitude defined as in equation 1 5, the differential of equation 1 5 may 
be expressed as 

g^ • dH = g • dZ. (16) 


This expression, introduced into equation 5, will reduce the number of variables prior to 
its integration, thereby leading to an expression for computing pressure as a function of 
geopotential height. 

The inverse-square law of gravitation provides an expression for g as a function of altitude 
with sufficient accuracy for most model-atmosphere computations: 
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( 17 ) 



where r Q is the effective radius of the earth at a specific latitude as given by Lambert’s 
equations (List, 1968). Such a value of r Q takes into account the centrifugal acceleration 
at the particular latitude. For this Standard, the value of r Q is taken as 6356.766 km, and 
is consistent with the adopted value of g Q = 9.80665 m/s 2 for the sea-level value of the 
acceleration of gravity. The variation of g as a function of geometric altitude is depicted 
in figure 2. 



Figure 2. Acceleration of gravity versus geometric altitude. 
Integration of equation 1 5, after substitution of equation 17 for g, yields 
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or 


r • H 

Z = — 09) 

r * r 0 - H 

where r = g 0 /gj ) = 1 m'/m. 

Differences between geopotential altitudes obtained from equation 18 for various values 
of Z, and those computed from the more complex relationship used in developing the 
U.S . Standard Atmosphere , 1962, (COESA, 1962) are small. For example, values of H 
computed from equation 18 are approximately 0.2, 0.4, and 33.3 m greater at 90, 120, 
and 700 km, respectively, than those obtained from the relationship used in the 1962 
Standard. In the 1975 Standard, geopotential altitude is used explicitly only at heights 
below 86 geometric kilometers. 

The transformation from Z to H (equation 16) in the development of the pressure-height 
relationship, for heights between the surface of the earth and 86-km altitude, makes it 
necessary to define the altitude variation of T and M in terms of H. It is convenient, there- 
fore, to determine the sea-level value of M, as well as the extent of any height dependence 
of this quantity, between the surface of the earth and 86-km altitude. Then, for this low- 
altitude regime, the two variables, T and M, are combined with the constant M 0 into a 
single variable T M , which is then defined as a function of H. 


Mean Molecular Weight 

The mean molecular weight, M, of a mixture of gases is by definition 

2 (n . • M .) 

M = (20) 


where n. and M. are the number density and defined molecular weight, respectively, of the 
i th gas species. In that part of the atmosphere, between the surface of the earth and about 
80-km altitude, mixing is dominant, and the effect of diffusion and photochemical proces- 
ses upon M is negligible. In this region, the fractional composition of each species is as- 
sumed to remain constant at the defined value, F. , and M remains constant at its sea-level 
value, M q . For these conditions, n. is equal to the product of F. times the total number 
density, N, so that equation 20 may be rewritten as 


2 [F j • N(Z) • MjJ ^(F i • Mj) 

M = M = : = 1 

0 2[F.*N(Z)] 2F j 


( 21 ) 


The right-hand element of this equation results from the process of factoring N(Z) out of 
each term of both the numerator and the denominator of the preceding fraction, so that, 
in spite of the altitude dependence of N, M is seen analytically to equal M 0 over the entire 
altitude region of complete mixing. When the defined values of F. and ML (from table 2) 
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are introduced into equation 21 , M Q is found to be 28.9644 kg/kmol. At 86 km (84.852 km’), 
however, the defined value of atomic-oxygen number density (8.6 X 10 16 m' 3 ) is seen (in 
Appendix A) to lead to a value of M = 28.9522 kg/kmol, about 0.04 percent less than M Q . 

To produce a smooth transition from this value of M to M 0 , the altitude profile of M has 
been arbitrarily defined at intervals of 0.5 km', for altitudes between 79.006 and 84.852 km', 
in terms of the ratio M/M Q as given in table 7. These ratio values have been interpolated 
from those initially selected for intervals of 0.5 geometric kilometers between 80 and 86 km 
to satisfy the boundary conditions of M = M 0 = 28.9644 at 80 km, and M = 28.9522 at 
86 km, and to satisfy a condition of smoothly-decreasing first differences in M within the 
height interval, 80 to 86 km. 


Table 7 

Molecular-weight Ratio versus Geopotential and Geometric Altitudes in Meters 


H 

Z 

M/M 0 

Z 

H 

M/M 0 

79000 

79994.1 

1.000000 

80000 

79005.7 

1.000000 

79500 

80506.9 

0.999996 

80500 

79493.3 

0.999996 

80000 

81019.6 

0.999988 

81000 

79980.8 

0.999989 

80500 

81532.5 

0.999969 

81500 

80468.2 

0.999971 

81000 

82045.4 

0.999938 

82000 

80955.7 

0.999941 

81500 

82558.6 

0.999904 

82500 

81443.0 

0.999909 

82000 

83071.5 

0.999864 

83000 

81930.2 

0.999870 

82500 

83584.8 

0.000822 

83500 

82417.3 

0.999829 

83000 

84098.0 

0.999778 

84000 

82904.4 

0.999786 

83500 

84611.4 

0.999731 

84500 

83391.4 

0.999741 

84000 

85124.8 

0.999681 

85000 

83878.4 

0.999694 

84500 

85638.4 

0.999679 

85500 

86000 

84365.2 

84852.0 

0.999641 

0.999578 


These arbitrarily-assigned values of M/M Q may be used for correcting a number of parame- 
ters of this Standard, if the tabulations are to correctly fit the model in the fifth, and per- 
haps in the fourth, significant figures within this height region. This after-the-fact correc- 
tion is required because these values of M/M 0 were not included in the program used for 
computing the tables of this Standard below 86 km, and hence, the tabulations of some of 
the properties may show a discontinuity of up to 0.04 percent between 85.5 and 86 km. 
This situation exists particularly for four properties in addition to molecular weight: 
kinetic temperature, total number density, mean-free path, and collision frequency. For 
these five parameters, the discrepancy in the detailed tables between 80 and 86 km can be 
readily remedied by a simple multiplication or division: Tabulated values of M, T, and L 
must be multiplied by the corresponding values of M/M 0 from table 7; tabulated values of 
N and v must be divided by the corresponding values of M/M 0 . 
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Three other properties-dynamic viscosity, kinematic viscosity, and thermal conductivity 
(which are tabulated only for heights below 86 km)-have similar discrepancies for heights 
immediately below 86 km. These values are not so simply corrected, however, because of 
the empirical nature of their respective defining functions. Rather, these quantities must 
be recalculated in terms of a suitably-corrected set of values of T, if the precisely correct 
values are desired for geometric altitudes between 80- and 86-km altitude. 


Molecular-scale Temperature versus Geopotential Altitude (0 to 84.3520 km') 

The molecular-scale temperature, T M (Minzner et al., 1958), at a point is defined as the 
product of the kinetic temperature, T, at that point times the ratio, M Q /M, where M is the 
mean molecular weight of air at that point and M Q is the sea-level value of M discussed 
above (see Appendix C). Analytically, 


T 


M 



( 22 ) 


(When T is expressed in the Kelvin scale, T M is also expressed in that scale.) 

The principle virtue of the parameter T M is that it combines the variable portion of M with 
the variable T into a single new variable, in a manner somewhat similar to the combining of 
the variable portion of g with Z to form the new variable H. When both of these transfor- 
mations are introduced into equation 5, and when T M is expressed as a linear function of 
H, the resulting differential equation has an exact integral. Under these conditions, the 
computation of P versus H becomes a simple process not requiring numerical integration. 
Traditionally, standard atmospheres have defined temperature as a linear function of height 
to eliminate the need for numerical integration in the computation of pressure versus 
height. This Standard follows the tradition to heights up to 86 km, and the function T M 
versus H is expressed as a series of seven successive linear equations. The general form of 
these linear equations is 

b + ^m, b * ~ H b ) (23) 


with the value of subscript b ranging from 0 to 6 in accordance with each of seven succes- 
sive layers. The value of T M for the first layer (b = 0) is 288. 1 5 K, identical to the sea- 
level value of T, since at this level M = M 0 . With this value of T M defined, and the set of 
six values of H b and the six corresponding values of L M b defined in table 3, the function 
T m of H is completely defined from the surface to 84.852 km' (86 km). A graph of this 
function compared to the similar function of the 1962 Standard is shown in figure 3. From 
the surface of the earth to the 51-km' altitude, this profile is identical to that of the 1962 
Standard. The profile from 51 to 84.852 km' was selected in accordance with present-day 
data, and abbreviated tables of thermodynamic properties of the atmosphere based upon 
this temperature-height profile were published by Kantor and Cole (1973). 
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MOLECULAR-SCALE TEMPERATURE, K 


Figure 3. Molecular-scale temperature versus geopotential 
altitude. 

Kinetic Temperature versus Geometric Altitude (0 to 1000 km) 

Between the surface of the earth and 86-km altitude, kinetic temperature is based upon the 
defined values of T M . In the lowest 80 km of this region, where M is constant at M Q , T is 
equal to T M in accordance with equation 22. Between 80 and 86 km, however, the ratio 
M/M 0 is assumed to decrease from 1 to 0.9995788, as indicated in table 7, such that the 
values of T correspondingly decrease from those of T M . Thus, at Z 7 = 86 km, a torm ot 
equation 22 shows that T ? has a value 186.8673 K, that is, 0.0787 K smaller than that of 
T M at that height. 

At heights above 86 km, values of T M are no longer defined, and geopotential is no longer 
the primary argument. Instead, the temperature-altitude profile is defined in terms of four 
successive functions, each of which is specified in such a way that the first derivative of T 
with respect to Z is continuous over the entire altitude region, 86 to 1000 km. These four 
functions begin successively at the first four base heights, Z b , listed in table 5, and are de- 
signed to represent the following conditions: 

1. An isothermal layer from 86 to 91 km; 

2. A layer in which T(Z) has the form of an ellipse from 91 to 110 km; 
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3. A constant, positive-gradient layer from 1 10 to 120 km; and 

4. A layer in which T increases exponentially toward an asymptote, as Z increases 
from 120 to 1000 km. 


86 to 91 km 

For the layer from Z ? = 86 km to Z g =91 km, the temperature-altitude function is defined 
to be isothermally linear with respect to geometric altitude, so that the gradient of T with 
respect to Z, is zero (see table 5). Thus, the standard form of the linear function, which is 

T = T b +L K,I. -< Z - Z b> <24) 

degenerates to 


T = T ? = 186.8673 K 


(25) 


and by definition 


— = 0.0 K/km 
dZ 


(26) 


The value of T 7 is derived from one version of equation 22 in which T M is replaced by 
T m 7 = 186.946 from equation 23 or from table 3, and M/M Q is replaced by the value 
0.9995788 from table 7. Thus, T ? = 186.8673 K. Since the kinetic temperature, T, is 
defined to be constant for the entire layer, Z ? to Z g , the temperature at Z g is T g = T ? 
= 1 86.8673 K, and the gradient, dT/dZ, at Z g is L K g = 0.0 K/km, the same as for L K 7 


91 to 1 10 km 

For the layer Z g = 91 km to Z g = 1 10 km, the temperature-altitude function is defined to 
be a segment of an ellipse expressed by 



(27) 


where 

T c = 263.1905 K, derived in Appendix B, 
A = -76.3232 K, derived in Appendix B, 

a = -19.9429 km, derived in Appendix B, 
and Z is limited to values from 91 to 110 km. 
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Equation 27 is derived in Appendix B from the basic equation for an ellipse, to meet the 
values of T g and L K derived above, as well as the defined values T 9 = 240 K and 
L k = 1 2 K/km, for Z 9 = 110 km. With these restraints, the values of T c , a, and A are 
found to be those cited above. 


The expression for dT/dZ related to equation 

dT -A / Z ~ Z 8 \ 

dZ a \ a / 

1 10 to 120 km 

For the layer Z 9 = 1 10 km to Z 10 = 120 km, T(Z) has the form of equation 24, where 
subscript b is 9, such that T b and L K b are, respectively, the defined quantities T 9 and L K 9 
(see Category III Constants and table 5 respectively), while Z is limited to the range 1 10 to 
1 20 km. Thus, 

T = T 9 + L k 9 (Z- Z 9 ) (29) 

and 



— = L„ 0 = 12.0 K/km (30) 

dZ K - 9 

Since dT/dZ is constant over the entire layer, L K 10 , the value of dT/dZ at Z 10 , is identical 
to L k 9 (that is, 1 2 K/km) while the value of T, 0 at Z, 0 is found from equation 29 to be 
360 K. 


120 to 1000 km 

For the layer Z 10 = 120 to Z 12 = 1000 km T(Z) is defined to have the exponential form 
(Walker, 1965) 

T = T.- (T_- T, 0 )- exp(-X* & (31) 

such that 


because 

and 


dT 

dZ 


^•CT.-T 10 ) 



• exp (- X •£) 


X = L K9 /(T.-T 10 ) = 0.01875, 

% = *(Z) = (Z-Z, 0 ) * (r Q +Z, 0 )/(r 0 + Z) 


(32) 
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In the above expressions, equals the defined value 1000 K. A graph of T versus Z from 
0- to 1000-km altitude is given in figure 4. The upper portion of this profile was selected 
to be consistent with satellite drag data (Jacchia, 1971), while the mid-portion, particularly 
between 86 and 200 km, and to some extent in the region 200 to 450 km, was selected* to 
be consistent with observed temperatures and satellite observations of composition data 
(Hedin et ah, 1972). 



Figure 4. Kinetic temperature versus geometric altitude. 

COMPUTATIONAL EQUATIONS 

The tables of this Standard have been computed in two height regions, 0 to 86 km 
(84.852 km'), and 86 to 1000 km, because the computations for each region are based 
on compatible but different sets of initial conditions. These two different sets of initial 
conditions lead to two different computational procedures. Consequently, the following 
discussion of computational equations, which is presented according to a series of atmo- 
spheric parameters, does not necessarily flow in the order in which the calculation is 
actually performed for each altitude region. The equations used for computing the various 


*Minzner, R. A., C. A. Reber, K. S. W. Champion, F. T. Huang, O. K. Moe, A. O. Nier, G. R. Swenson, S. P. Zimmerman, 
“The 1975 Standard Atmosphere Above 86-km Altitude: Recommendations of Task Group II to COESA,” 1974, to be 
published as NASA SP. 
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properties of the atmosphere for altitudes below 86 km are, with certain noted exceptions, 
equivalent to those used in the 1962 Standard, and the various equations involving T M stem 
from expressions used in the ARDC Model Atmosphere, 1956 (Minzner, 1956). 


Pressure 


Three different equations are used to compute pressure, P, in various height regimes of this 
Standard. One of these equations applies to heights above 86 km, while the other two apply 
to the height regime from the surface of the earth up to 86-km altitude, within which re- 
gime the argument of the computation is geopotential. Consequently, expressions for com- 
puting pressure as a function of geopotential altitude stem from the integration of equation 
5 after replacing g • dZ by its equivalent g^ dH from equation 16, and after replacing the 
ratio M/T by its equivalent M 0 /T 0 in accordance with equation 22. Two forms result from 
this integration— one is for the case when L M for a particular layer is equal to zero, and the 
other when the value of L M b is not zero. The latter of these two expressions is 


P = Pu 


V M o 


M, b 


R* * L 


M, b 


T M,b + L M,b H b ^ 


and the former is 


(33a) 


P = 


P, • exp 


g; • M 0 (H- H b ) 


R* • T. 


M, b 


(33b) 


Equation 33a is used for layers associated with values of subscript b equal to 0, 2, 3, 5 and 
6; equation 33b is used for layers associated with values of subscript b equal to 1 and 4. 

In these equations, g^, M Q , and R* are each defined single-valued constants, while L M b and 
EL are each defined multi-valued constants in accordance with the value of b as indicated 

b 

in table 3. In each equation, H may have values ranging from H b to H b+1 . The quantity 
T m is a multi-valued constant listed in table 3 with values derived from equation 23 in 
accordance with the several values of b and the corresponding defined values of L M b and 
FL . The reference-level value of P u for b = 0 is the defined sea-level value, P^ = 101 .3250 

b b ’0 

kPa (equivalently 101325 N/m 2 or 1013.25 mbar). The values of P b for b = 1 through 
b - 6 are obtained from the application of the appropriate member of the pair of equations 
33a and 33b for the case when H = H b+1 . 

These two equations applied successively yield the pressure for any desired geopotential 
altitude from sea level to H ? , where H 7 is the geopotential altitude corresponding to the 
geometric altitude Z ? = 86 km. Pressures for H from 0 to -5 km' may also be computed 
from equation 33a when b = 0. 
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For Z equal to 86 km and above the value of pressure is computed as a function of geome- 
tric altitude, Z, and involves the altitude profile of kinetic temperature, T, rather than that 
of T m , in an expression in which the total pressure, P, is equal to the sum of the partial 
pressures for the individual species as expressed by equation 3. Thus, for Z = 86 to 1000 km, 

P = 2P. = In. • k • T = — 5 (33c) 

1 ’ N. 

A 


In this expression, 

k = the Boltzmann constant, defined in Category I, 

T = T(Z) defined in equations 25, 27, 29, and 31 for successive layers, 

2m = the sum of the number densities of the individual gas species comprising the 
atmosphere at altitude Z above 86 km, as described below. 

Neither m, the number densities of individual species, nor 2m, the sum of the individual 
number densities, is known directly. Consequently, pressures above 86 km cannot be com- 
puted without first determining m for each of the significant gas species. 

Number Density of Individual Species 

In the height region of complete mixing (0 to 80 km), n., the number density of any parti- 
cular major gas species varies with altitude in accordance with the altitude variation of the 
total number density, N. For this height region, the value of the i th species is, therefore, 
given by 

n i = Fj • N (34) 

where F. is the constant fractional volume coefficient given for each species in table 2. 

Since the values of N listed in the detailed tables of this Standard are not completely con- 
sistent with the basic definition of the model between 80 and 85 km, as previously dis- 
cussed, values of n. calculated from tabulated values of N for this limited height regime, 
must be corrected by dividing by the appropriate values of M/M Q from table 7. At altitudes 
above 86 km, equation 34 no longer applies, since the model assumes the existence of vari- 
ous processes which lead to particular differing height variations in the number-density 
values of several individual species, N 2 , O, 0 2 , Ar, He, and H, each governed by equation 6. 
Ideally, the set of equation 6, each member of which is associated with a particular species, 
should be solved simultaneously, since the number densities of all the species are coupled 
through the expressions for molecular diffusion which are included in equation 6. Such a 
solution would require an inordinate amount of computation, however, and a simpler 
approach was desired. This was achieved with negligible loss of validity by some simplifying 
approximations, and by calculating the number densities of individual species one at a time 
in the order n(N 2 ), n(O), n(0 2 ), n(Ar), n(He), and n(H). For all species except hydrogen 
(which is discussed later), we divide equation 6 by n. and integrate directly, to obtain the 
following set of simultaneous equations, one for each gas species: 
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n i = "i, 7 


• exp 


f(Z) + 


D. + K 
1 


In this set of equations, 

n = the set of species-dependent, number-density values for Z = Z ? = 86 km, 

one member for each of the five designated species, as derived in Appen- 
dix A and listed here, 

N 1 . 1 29794 X 1 0 20 m' 3 

O 8.6 X 10 16 nr 3 

0 2 3.030898 X 10 1 9 nr 3 

Ar 1.351400 X 10 18 m 3 

He 7.5817 X 10 10 nr 3 

T 7 = 1 86.8673 K, the value of T at Z ? , as specified in equation 25, 

T = T(Z) defined by equations 25, 27, 29, and 31 for the appropriate altitude 

regions, 

f(Z) = the function written as equation 36 below, 

v./(D. + K) = the set of empirical functions written as equation 37 below 
For f(Z) we have 


D. + K 

i 


M • K a i R 
M. + + 

1 D g 


where 


dT/dZ 


= D.(Z) as defined by equation 8 for the ith species, 

= K(Z) as defined by equations 7a, 7b, and 7c, 

= the molecular weight of the ith species as defined in table 2, 

= the thermal diffusion coefficient for the ith species as defined in table 4, 
= one of equations 26, 28, 30, or 32, as appropriate to the altitude region, 
= M(Z), with special considerations mentioned below. 


For [v./(D.+ K)] we have the following set of empirical expressions. 


D + K 


= Qj • (Z - Uj) 2 • exp [-Wj • (Z- U,) 3 ] 


+ q ; • (u ; - Z) 2 • exp [-Wj • (u ; - Z) 3 ] 
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The set of expressions represented by equation 37, while representing a function of both D. 
and K, does not directly involve calculated values of either of these coefficients. Rather, it 
involves a series of six other coefficients which, for each of four species, have been empiri- 
cally selected to adjust the number-density profile of the related species to particular values 
in agreement with observations. The defined values of the six sets of species-dependent 
coefficients— q. , Q., u., U. , w., and W. used in equation 37— are listed in table 6. The values 
of q. and U. were selected so that for 0 2 , Ar, and He, the quantity v./(D. + K) becomes 
zero at exactly 86 km. For atomic oxygen, however, all six of these coefficients contribute 
to maximizing this quantity for Z = 86 km. 


Molecular Nitrogen 

Molecular nitrogen (N 2 ) is the first species for which n is calculated. On the average, the 
distribution of N 2 is close to that for static equilibrium, and hence, for this species, we 
may neglect the transport velocity, thereby eliminating the term [v./(D. + K)] from that 
version of equation 35 applying to N 2 . This species is dominant above and below the 
turbopause, and its molecular weight is close to the mean molecular weight in the lower 
thermosphere, where mixing still dominates the distribution process. The effect of mixing 
up to 100-km altitude is approximated therefore by two additional adjustments to equation 
35 as applied to N 2 . Both adjustments are implicit in f(Z); these are neglecting K and re- 
placing M. by the mean molecular weight M which, for the altitude region 86 to 100 km,*, 
is approximated by M Q . With these three adjustments, that version of equation 35 applying 
to N 2 reduces to 


n(N 2 ) 


n(N 7 ) 7 * exp 




(38) 


where 

M = M Q for Z < 100 km, and 
M = M(N 2 ) for Z> 100 km. 

Figure 5 shows a graph of n(N 2 ) versus Z. 

Species O, 0 2 , A r, and He 

As noted above, after the calculation of n(N 2 ) has been performed, the values of n. for the 
next four species are calculated from equation 35 in the order O, 0 2 , Ar, and He. In the 
case of O and 0 2 , the problem of mutual diffusion is simplified by considering N 2 as the 
stationary background gas (as described in the previous section). For Ar and He, which 
are minor constituents in the lower thermosphere, it is more realistic to use the sum of 
the number densites of N 2 , O, and 0 2 as the background gas in evaluating the molecular- 
diffusion coefficient, D., and the mean-molecular weight, M, except below 100 km where 


26 




Figure 5. Number density of individual species and total 
number density versus geometric altitude. 


M is taken to be the sea-level value, M 0 . This latter choice is to maintain consistency with 
the method for calculating n(N 2 ). 

In equation 37, defining [v./(D. + K)] , the coefficients q ; , u., U., w., and W., which 

(except for q.) are constant for a particular species, are each adjusted such that appropriate 
densities are obtained at 450 km for O and He, and at 1 50 km for O, 0 2 , He, and Ar. The 
constant q., and hence the second term of equation 37, is zero for all species except atomic 
oxygen, and is also zero for atomic oxygen above 97 km; the extra term for atomic oxygen 
is needed below 97 km to generate a maximum in the density-height profile at the selected 
height of 97 km. This maximum results from the increased loss of atomic oxygen by re- 
combination at lower altitudes. The flux terms for O and 0 2 are based on, and lead (quali- 
tatively) to the same results as those derived from the much more detailed calculations by 
Colegrove et al. (1965) and Keneshea and Zimmerman (1970) and discussed in Appendix 
D. 

A further computational simplification is realized above 1 15 km where the eddy-diffusion 
coefficient becomes zero. For these altitudes, the set of expressions represented by equa- 
tion 35 becomes uncoupled, and each member reduces to a form where the integration is 
performed only on the sum of three terms: 

1 . The barometric term for the particular species (that is, the right-hand side of 
equation 5), 
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2. The thermal-diffusion term (cx./T) * (dT/dZ), 

3. A simplified velocity term. 

In the case of O, 0 2 , and Ar, the thermal-diffusion term is zero. Also, as may be shown, 
the velocity term, [ v. /(D. + K)] , becomes small above 120 km and, with the exception of 
atomic hydrogen, each species considered is nearly in diffusive equilibrium at these heights. 
For the present model, however, this situation becomes exactly true only at altitudes above 
150 km. 

The altitude profile of number density for each of the species O, 0 2 , Ar, and He is given 
in figure 5, along with that for N 2 . 


Atomic Hydrogen 


For various reasons, the height distribution of the number density of atomic hydrogen, 
n(H), is defined only for heights from 1 50 to 1000 km. Below 150 km, the concentration 
of H is negligible compared with the concentrations of O, 0 2 , Ar, and He. The defining 
expression for n(H), like the expression for n(N 2 ), n(O), and so on, is derived from equa- 
tion 6. The solution for n(H), however, is expressed in terms of the vertical flux, n(H) 

• v(H) represented by 0, rather than in terms of v(H), because it is the flux which is con- 
sidered known for H. In this model, only that contribution to 0 due to planetary escape 
from the exosphere is considered. 


Since K is zero for the altitude region of interest, the particular version of equation 6 ap- 
plied to H is correspondingly simplified, and one possible solution to the resulting expres- 
sion is 



pz <t> / 1 

p \l+a(H) 

n(H) = 

n(H)|1 X °( H ) ( T 

L ii y 

— J • (exp r) * dZ 

n / 


1 + a(H) 


(exp -r) 


(39) 


where 


n(H) n 

D(H) 


<P 


T 

T 

li 


8.0 X 10 10 m' 3 , the number density of H at Z n = 500 km, as defined in 
Category III Constants, 

the molecular diffusion coefficient for hydrogen given by equation 8 in which 
the values of a. and b { are as defined in table 4, 

7.2 X 10 11 m' 2 • s' 1 , the vertical flux of H, as defined in Category III Constants, 
T(Z) as specified by equation 31, 

999.2356 K, the temperature derived from equation 31 for Z = Z n , 
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a(H) = the thermal diffusion coefficient for H, -0.25 (dimensionless), as defined in 
table 4, 

t = t( Z) defined in equation 40 


T 



g • M(H) 

R* • T 


dZ 


( 40 ) 


Because D(H) becomes very large compared with 0 for heights above 500 km, the value of 
the integral term in equation 39 can be neglected at these heights, and atomic hydrogen is 
then essentially in diffusive equilibrium. Figure 5 depicts the graph of n(H) as a function 
of Z, along with those of other species. 


Mean Molecular Weight (Above 86 km) 

Equations 35 through 39 permit the calculation of the number densities of the species N 2 , 
O, 0 2 , Ar, He, and H for heights above 1 50 km, and of the first five of these species for 
heights between 86 and 150 km, where n(H) is insignificant compared with n(N 2 ). These 
number densities permit the calculation of several atmospheric parameters in the height 
region 86 to 1000 km. The first is mean molecular weight using equation 20. These values 
of M, along with those implicit in table 7 for Z from 80 to 86 km, plus the invariant value, 
M Q , for heights from 0 to 80 km, are shown in figure 6. 



Figure 6. Mean molecular weight versus geometric altitude. 
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Pressure and Mass Density (Above 86 km) 

The number densities of the several species also permit us now to compute total pressure 
and mass density for heights from 86 to 1000 km, using equation 33c, and one version of 
equation 1. Figure 7 depicts these values as well as those for heights below 86 km computed 
from equations 33a and 33b. Finally, these individual number densities permit the calcula- 
tion of total number density, N = Zn., at heights from 86 to 1000 km. 

DENSITY, kg/m 3 

10" 14 10 ' 12 10" 10 10 8 10 ' 6 10 4 10" 2 10 ° 



Figure 7. Total pressure and mass density versus 
geometric altitude. 


Total Number Density 


From equations 2, 22, and 33c, it is apparent that total number density, N, the number of 
neutral gas particles per unit volume of the atmosphere, may be expressed in any one of 
the following three equivalent forms: 


M o • n a • p _ n a - p 
T M • R* • M R* • T 


(41) 
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The three forms are selected to satisfy three types of calculations: 

1. Those depending upon values of T M , 

2. Those depending upon values of T, and 

3. Those depending upon values of 2n.. 

This format will be followed in specifying the computational equations, insofar as possible, 
for the balance of the quantities discussed in this document. A graph of the altitude varia- 
tion of total number density is presented in figure 5, along with the number densities of 
individual species. That portion of equation 41 involving T M is of particular interest in 
calculating N for heights from 0 to 86 km. 


Mass Density 

From equations 1, 22, and 33c, one may write the following three forms of computational 
equations for mass density, p: 

p • M 0 p ■ M 2 ("j ‘ M i) 

P = = = (42 

R * • T M R* • T N a 

The altitude-dependent variations of this quantity are shown in figure 7 along with those 
of pressure. 


Mole Volume 

Mole volume, v , of air is defined as the volume of one mole of air, where one mole of air 
is the amount consisting of a number of neutral particles equal to N A . In S.I. units, the 
quantity v m should specify the number of cubic meters containing 1 kmol of air. Since M 
has the dimensions of kg/kmol, and p has the dimensions of kg/m 3 , the ratio M/p, with the 
units m 3 /kmol, provides the definition of mole volume. In terms of equations 1 , 22, and 
33c, this ratio may be equated to the following series of expressions: 


v 


rn 


R* • M • T 


M 


M 


0 


P 


R* • T _ N A 
P Sn; 


(43) 


This quantity, while not tabulated in this Standard, is shown graphically in figure 8. 


Scale Height 

Pressure Scale Height 

The quantity R* • T/(g • M), which has dimensions of length, is commonly associated with 
the concept of scale height and is the defining form of pressure scale height, H p , used in 
this model. Thus, equations 20 and 22 may be written: 
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Figure 8. Mole volume versus geometric altitude. 

R* • T m R* • T R* • T • £ nj 

IT = = = (44 

p g • M q g • M g • Sfm • M ; ) 

The reciprocal of this quantity, which appears on the right-hand side of equation 5, is seen 
to equal the slope of the function £nP versus Z at height Z in the regions where hydrostatic 
equilibrium or diffusive equilibrium holds. In the present model, this condition is true for 
heights below 80 km (complete mixing), and is essentially true above approximately 
120 km where diffusive equilibrium is nearly satisfied, and where each individual species 
is governed by equation 4. 

In the 80- to 120-km region, where the transition from a completely mixed atmosphere to 
one in diffusive equilibrium takes place, the situation is complicated by the competition 
between three processes-molecular diffusion, eddy diffusion, and dissociation of molecu- 
lar oxygen. These processes result in a vertical transport, such that equations 4 and 5 are 
no longer exactly true in this 40-km layer. Since molecular nitrogen is the dominant spec- 
ies in this altitude range, however, and since this species has a zero transport velocity in 
this model, the pressure scale height is still a good indicator of the rate of change of pres- 
sure in this height region. 
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It should be noted that equations 4 and 5 also become invalid at very high altitudes (the 
exospheric region) due to the infrequent collisions between neutral particles. Thus, in this 
region, the significance of H p as a measure of dCnP/dZ again loses validity. 

In equation 44, both g and T M or all three of g, T, and M are functions of Z such that H p 
is the local value of geometric pressure scale height. This quantity, which is the particular 
scale height tabulated in this Standard, and which is plotted in figure 9, is frequently but 
incorrectly associated with the altitude increment over which the pressure decreases by 
exactly a factor of 1/e. The conditions necessary for the pressure to decrease by exactly 
that factor over an altitude increment of a single pressure scale height would be for the 
variabies T, g, and M all to remain constant over that altitude interval. Since g may never 
be constant over any altitude interval, this particular concept of pressure decrease can 
rarely, if ever, apply exactly to H p . 



Figure 9. Pressure scale height and density scale 
height versus geometric altitude. 


With the replacement of g • dZ in equation 5 by g Q • dH in accordance with equation 16, 
and with the replacement of T/M in equation 5 by T M /M 0 in accordance with equation 22, 
as in the development of equations 33a and 33b, an expression for the slope of the func- 
tion CnP versus geopotential height is developed. The negative reciprocal of the expression 
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for this slope is called geopotential pressure scale height, H p , and it has the same dimensions 
as geopotential height. The quantity Hp is implicit in the computational pressure equations 
33a and 33b. This concept of scale height may be expressed as follows: 


/ -dKnP V 1 R * ' T M R* • T 
\ dH / " gy • M 0 ‘ gy • M 


(45) 


Within a layer of constant T M , or of constant ratio T/M, the pressure does, in fact, decrease 
by exactly 1/e of its initial value when the geopotential altitude increases by exactly one 
geopotential scale height. This is precisely the situation expressed by equation 33b which, 
in terms of Hp, may be rewritten as 


P = 



(46) 


Equation 45 is rigorously valid only below 80 km, while equation 46 applies only within 
isothermal layers below 80 km. 


Density Scale Height 


Because of the relationship between H p and the slope of CnP versus Z it is convenient to 
apply the name geometric density scale height, H p , to the negative reciprocal of the slope 
of Cnp versus Z. Using the equation of state (equation 1) to relate dCnp to dCnP, one may 
define 


H 


p 


H 


P 


1 + H p 



1 + H p 



dCnM 

dZ 


(47) 


The relationships implied between H p and dCnp/dZ are subject to the same limitations as 
those between H p and dCnP/dZ expressed above, that is, H p is only an approximation to 
(dfinp/dZ)" 1 between 80 and 120 km and in the exosphere, where the approximation be- 
comes increasingly invalid with increasing altitude. 

Within these limitations, it is apparent that in layers where T M does not change with 
changing altitude, that is, where (d£nT M /dZ) = 0, H p is equal to H p . Within such layers, 
the slope of finP versus Z at any particular altitude Z is identical to the slope of CnP versus 
Z, but neither P nor p decreases by a factor of 1 /e when altitude increases by an amount 
H P = H p . 

While density scale height is not tabulated in this Standard, values of this quantity are 
shown graphically with those of H p in figure 9. 
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If, however, a geopotential density scale height, FT, is defined as the negative reciprocal of 
the slope of fnp with respect to geopotential, one may write 




H p 


t / d^nT dCnM \ 

Hp y dH dH J 


1 + 


(48) 


This quantity is equal to Hj, in a layer within which T M is constant, and for this condition, 
IT specifies the altitude change for which p changes by a factor of 1/e. The quantity FP 
permits the writing of a density-altitude equation analagous to the pressure-altitude equa- 
tion, equation 33b, but in terms of FP : 


P 



(49) 


As in the case of equation 45, equation 49 is rigorously valid only below 80 km, while 
equation 49 applies only within isothermal layers below 80 km. 


Mean Air-particle Speed 

The mean air-particle speed, V, is the arithmetic average of the speeds of all air particles in 
the volume element being considered. All particles are considered to be neutral. For a 
valid average to occur, there must be a sufficient number of particles involved to represent 
mean conditions. Pressure and temperature gradients within the volume must also be neg- 
ligible. The analytical expression for V is closely related to that for the speed of sound, 
and is proportional to the ratio T/M. Thus, in terms of equations 20 and 22, 


8 • R* • T m 

1/2 

8 * R* * T 

1/2 

8 • R* • T • Zn [ 

7T • M 0 


71 * M 


7r ■ 2(ni • M ; ) 


The variation of particle speed with geometric altitude is shown in figure 10. 


Mean Free Path 


The mean free path, L, is the mean value of the distances traveled by each of the neutral 
particles, in a selected volume, between successive collisions with other particles in that 
volume. As in the case of V, a meaningful average requires that the selected volume be 
big enough to contain a large number of particles. The computational form for L is 


L = 


2 1/2 • R* • M • T M 
2tt • N a • a 2 • M 0 • P 


2 1/2 . r* . j 
2n • N a • a 2 • P 



(51) 
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MEAN PARTICLE SPEED, m/s 

Figure 10. Mean air-particle speed versus geometric 
altitude. 

where a is the effective collision diameter of the mean air molecules. The adopted value of 
a, that is, 3.65 X 10' 10 m is suitable for that part of the atmosphere below about 86 km, 
which is dominated by N 2 and 0 2 . Above this height, the value of o, which depends upon 
composition in a complicated manner, begins to change significantly so that tabulations 
with four significant figures are no longer valid. At great altitudes, this expression for L 
is valid only under assumptions that hold M, T M , P, and a constant throughout the volume 
used. Figure 1 1 depicts the mean free path in terms of altitude. 


Mean Collision Frequency 

The mean collision frequency, v , is the average speed of the air particles within a selected 
volume divided by the mean free path, L, of the particles within that volume. That is, 


v 



(52) 


and in computation form: 
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MEAN FREE PATH, m 

Figure 1 1 . Mean free path versus geometric altitude. 


v = 4N a * o 
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A 


7 r • r 
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(53) 


= 4N^ • a 1 • 


n ■ P z * 2n- 


R* • T • S(n j • M ; ) 


1/2 


Note that a is again involved in this quantity, and hence v has limitations similar to those 
of mean free path. The foregoing expressions are taken to apply to neutral particles only, 
since no considerations involving charged particles are introduced for purposes of develop- 
ing the tables and graphs of this Standard. 

Figure 1 2 graphically displays the variation of collision frequency with altitude. See Mean 
Air-particle Speed of this document for a discussion of the assumptions under which equa- 
tion 53 is valid at great altitudes. 


Speed of Sound 

The expression adopted for the speed of sound, C g , is 



(54) 
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Figure 12. Mean collision frequency versus geometric altitude. 

where 7 is the ratio of specific heat of air at constant pressure to that at constant volume, 
and is taken to be 1 A exact (dimensionless), as defined in Category II Constants. Equa- 
tion 54 for speed of sound applies only when the sound wave is a small perturbation on 
the ambient condition. Calculated values of C g have been found to vary slightly from 
experimentally determined values. 

The limitations of the concept of speed of sound due to extreme attenuation are also of 
concern. The attenuation, which exists at sea level for high frequencies, applies to succes- 
sively-lower frequencies as atmospheric pressure decreases, or as the mean free path in- 
creases. For this reason, the concept of speed of sound (except for frequencies approaching 
zero) progressively loses its range of applicability at high altitudes. Hence, the listing of the 
values for speed of sound are not given for heights above 86 km. Figure 13 shows the varia- 
tion with altitude of the computed speed of sound. 

Dynamic Viscosity 

The coefficient of dynamic viscosity, p, is defined as a coefficient of internal friction 
developed where gas regions move adjacent to each other at different velocities. The 
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SPEED OF SOUND, m/s 

Figure 1 3. Speed of sound versus geometric altitude. 


following expression, basically from kinetic theory but with constants derived from ex- 
periment, is used for computation of this quantity: 




0 • t 3/2 

T + S 


(55) 


In this equation, 0 is an exact constant equal to 1.458 X 10 -6 kg/(s • m • K 1/2 ), and S is 
Sutherland’s constant, equal to 1 10.4 K (exact), both defined in Category II Constants. 
Because of the empirical nature of this equation, no attempt has been made to transform 
it into one involving T M . 

Equation 55 fails for conditions of very high and very low temperatures, and under condi- 
tions occurring at great altitudes. Consequently, tabular entries for coefficient of dynamic 
viscosity are not given for heights above 86 km. For these reasons, caution is also necessary 
in making measurements involving probes and other objects which are small with respect to 
the mean free path of molecules, particularly in the region of 32 to 86 km. 


The variation of dynamic viscosity with altitude is shown in figure 14. 
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Figure 14. Dynamic viscosity versus geometric altitude. 


Kinematic Viscosity 

Kinematic viscosity, 17, is defined as the ratio of the dynamic viscosity of a gas to the densi- 
ty of that gas, that is. 


V 



(56) 


Limitations of this equation are comparable to those associated with dynamic viscosity, and 
consequently, tabular entries of kinematic viscosity are also not given for heights above 
86 km. Figure 15 is a graphical representation of the variation of kinematic viscosity with 
altitude. 


Coefficient of Thermal Conductivity 

The empirical expression adopted for purposes of developing tabular values of the coeffi- 
cient of thermal conductivity, k ( , for heights up to the 86-km level is as follows: 

2.64638 X I O' 3 • T 3/2 

k, (5 

T + 245.4 X 10' (12/T) 
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Figure 1 5. Kinematic viscosity versus geometric altitude. 

This expression differs from that used in the U.S. Standard Atmosphere, 1962 (COESA, 
1962), in that the numerical constant has been adjusted to accommodate a conversion of 
the related energy unit from the temperature-dependent kilogram calorie to the invariant 
joule. Thus, the values of k t in units of J/(s • m • K) or W/(m • K) are greater than the 
values of k ( in units of kcal/(s • m • K) by a factor of exactly 4.1858 X 10 3 , when the 
kilocalorie is assumed to be the one for 1 5°C. Kinetic-theory determinations of thermal 
conductivity of some monatomic gases agree well with observation. For these gases, ther- 
mal conductivity is directly proportional to the dynamic viscosity. Modification of the 
simple theory has accounted in part for differences introduced by polyatomic molecules 
and by mixtures of gases. Tabular entry of values for coefficient of thermal conductivity 
is terminated at 86 km. The variation with height of this quantity is shown in figure 16. 


41 





COEFFICIENT OF THERMAL CONDUCTIVITY, W/(m . K) 
Figure 16. Thermal conductivity versus geometric altitude. 


SUMMARY OF TABULAR VALUES OF ATMOSPHERIC PROPERTIES 
Sea-level Values 

The sea-level values of 1 5 of the atmospheric properties discussed in this Standard are as 
follows: 


Symbol 


C s,0 

3.40294 

&0 

9.80665 

H P,0 

8.4345 

^t,0 

2.5326 

L o 

6.6328 

V m ,0 

2.3643 

M o 

2.89644 

N o 

2.5470 

P o 

1.01325 

T o 

2.8815 

v 0 

4.5894 

% 

1.4607 

^0 

1.7894 

*0 

6.9193 

Po 

1.2250 


Sea-level Value 

X 10 2 m/s 
m/s 2 
X 10 3 m 

X 10‘ 3 J/(s • m • K) or W/(m • K) 
X 10 8 m 
X 10 l m 3 /kmol 
X 10* kg/kmol 
X 10 2S m' 3 
X 10 s N/m 2 
X 10 2 K 
X 10 2 m/s 
X 10 s m 2 /s 
X 10 s kg/(m • s) 

X 10 9 s 1 

kg/m 3 
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The sea-level values for g, P, and T are defined quantities; the remainder are quantities cal- 
culated from the preceding equations. 

Conversion of Metric to English Units 

For those who have need to work in the English System of Units, the conversion factors 
listed in table 8 are applicable to the atmospheric parameters tabulated or shown graphi- 
cally in this Standard. For other transformations, see Mechtly (1973). 


Table 8 

Metric to English Conversion Factors for Properties of the 
U.S. Standard Atmosphere, 1975 (COESA, 1975) 


Symbol 

To Convert from 
Metric Units 

To English 
Units 

Divide by 

c s 

m/s 

ft/s 

3.048 X 10 1 * 

g 

m/s 2 

ft/s 2 

3.048 X 10 1 * 

H p 

m 

ft 

3.048 X 10 1 * 

k t 

W/(m • K) 

BTU/(ft • s • °R) 

6.226 477 504 X 10 3 

L 

m 

ft 

3.048 X 10 1 * 

V m 

m 3 /kmol 

ft 3 /lbmol 

6.242 796 057 X 10' 2 

M 

kg/kmol 

lb/lbmol 

1.000* 

N 

_3 

m 

ft' 3 

3.531 466 672 X 10 2 

P 

mbar 

in Hg (32°F) 

3.386 389 X 10 1 

T or T m 

K 

°R 

5/9* 

V 

m/s 

ft/s 

3.048 X 10 1 * 

n 

m 2 /s 

ft 2 / s 

9.290 304 X 10" 2 * 

M 

N • s/m 2 

lb/(ft • s) 

1.488 163 944 

V 

s' 1 

s' 1 

1.000 

9 

kg/m 3 

lb/ft 3 

1.601 846 3 X 10 1 


*Exact definition. 


Tables of Atmospheric Properties 

Detailed tables of the height-dependent values of the various atmospheric properties de- 
fined by this standard have been prepared elsewhere (COESA, 1975). Abbreviated versions 
of these tables are presented at the end of this section in tables 9 through 1 5. The first four 
of these tables, 9 through 12, list atmospheric properties as a function of the boundary 
heights of the seven layers adopted to define the temperature-height profile between the 
surface of the earth and 86-km altitude. These boundary heights are specified in integer 
multiples of one geopotential kilometer, although the equivalent values of geometric height 
are also given. 
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Tables 13, 14, and 15 list the atmospheric properties of this Standard for 14 selected alti- 
tudes in the 86- to 1000-km interval. These altitudes include the boundary heights of the 
four layers adopted to define the temperature-height profile in this height region. Because 
of the large interval, 120 to 1000 km, comprising the highest layer of this model, the values 
of the various atmospheric properties in tables 13, 14, and 15 have been given for nine addi- 
tional heights within this layer. 

Only three tables (13, 14, and 15) are used to list the properties of this Standard for heights 
from 86 to 1000 km, as compared with four tables for the height region 0 to 86 km. This 
is due to the fact that, in this Standard, none of the four quantities-speed of sound, dyna- 
mic viscosity, kinematic viscosity, and thermal conductivity -is defined for heights above 
86 km. Consequently, table 12, which lists the values of these properties below 86 km, 
has no counterpart for heights above 86 km. Table 1 1, listing number densities of five 
atmospheric gas species for heights from 0 to 86 km, has no counterpart in the detailed 
tables of the COESA document (COESA, 1975). 
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Table 9 

Temperature, Pressure, and Density versus Geopotential Altitude, 0 to 86 km 


Altitude 

Temperature 

Pressure 

Density 

H 

Z 

T 

t 

T m 

P 

P 

P/P 0 

P 

pIp 0 

km' 

km 

K 

°C 

K 

mbar 

tori 


kg/m 3 

- 

0,0000 

0.0000 

288.150 

15.000 

288 150 

1.013250 +3 

7,60000 +2 

1 .00000 +0 

! .224999 +0 

1.0000+0 

11.0000 

11.0190 

216.650 

-56.500 

216.650 

2.263206 +2 

1.69754 +2 

2.23361 -1 

3.639178 -1 

2.9708 - 1 

20.0000 

20.0631 

216.650 

-56.500 

216.650 

5.474889 +1 

4 10650 +1 

5.40330-2 

8.803480-2 

7.1865-2 

32.0000 

32.1619 

228.650 

-44.500 

228.650 

8.680187 +0 

6.51068 +0 

8.56667-3 

1 .322500-2 

1.0796-2 

47,0000 

47.3500 

270.650 

- 2.500 

270.650 

1.109063 +0 

8.31866-1 

1.09456-3 

1.427532-3 

1.1653-3 

51.0000 

51.4124 

270.650 

- 2.500 

270.650 

6.693887-1 

5.02083 -1 

6.60635 -4 

8.616049 -4 

7.0335-4 

71.0000 

71.8019 

214.650 

-58.500 

214.650 

3.956420-2 

2.96756-2 

3.90468-5 

6 421099 -5 

5.2417-5 

84.8520 

86.0000 

186,867 

-86.283 

186.946 

3.733836-3 

2.80061 -3 

3.68501 -6 

6.957879-6 

5.6799-6 
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Table 10 

Seven Secondary Atmospheric Properties versus Geopotential Altitude, 0 to 86 km 


Altitude 

Acceleration 
Due to 
Gravity 

Pressure 

Scale 

Height 

Number 

Density 

Mean 

Particle 

Speed 

Mean 

Collision 

Frequency 

Mean Free 
Path 

Mean 
Molecular 
Weigh l 

H 

Z 

8 

H P 

N 

V 

V 

L 

M 

km f 

km 

m/s 2 

km 

m- 3 

m/s 

s- 1 

m 

kg/ k mol 

0.0000 

0.0000 

9.80665 

8.4345 

2.5471 +25 

458.94 

6.9193 +9 

6.6328-8 

28.9644 

11.0000 

11.0190 

9.77274 

6.3636 

7.5669 +24 

397.95 

1.7824 +9 

2.2327 -7 

28.9644 

20.0000 

20.0631 

9.74504 

6.3817 

1.8305 +24 

397.95 

4.31 17 +8 

9.2295 -7 

28.9644 

32.0000 

32.1619 

9.70817 

6.7608 

2.7499 +23 

408.82 

6.6542 +7 

6.1438 -6 

28.9644 

47.0000 

47.3500 

9.66217 

8.0407 

2.9683 +22 

444,79 

7.8146 +6 

5.6918-5 

28.9644 

51.0000 

51.4124 

9.64992 

8.0509 

1.7915 +22 

444.79 

4,7166+6 

9.4303-5 

28.9644 

71.0000 

71.8019 

9.58881 

6.4258 

1.3351 +21 

396.11 

3.1303 +5 

1.2654-3 

28.9644 

84.8520 

86.0000 

9.54659 

5.6212 

1 .4473 +20 

396.67 

3.1667 +4 

1.1674 -2 

28.9522 


Table 1 1 

Number Densities of Five Species versus Geopotential Altitude, 0 to 86 km 


Altitude 

Number Densities 

11 

Z 

n(N 2 ) 

u(0) 

n(0 2 ) 

n(A) 

n(He) 

km 

km 

m- 3 

m 3 

in' 3 

nf 3 

m' 3 

0.0000 

0.0000 

1 .9888 +25 

- 

5.3353 +24 

2,3789 +23 

1.3346+20 

11.0000 

11.0190 

5.9082 +24 

- 

1 .5850 +24 

7.0671 +22 

3.9648+19 

20.0000 

20.0631 

1.4292 +24 

- 

3.8342 +23 

1.7096 +22 

9.5912 +18 

32.0000 

32.1619 

2.1471 +23 

- 

5.7599 +22 

2.5682 +21 

1.4408+18 

47.0000 

47.3500 

2.3176 +22 

- 

6.2174+21 

2.7722 +20 

1.5553+17 

51.0000 

51.4124 

1 .3988 +22 

- 

3.7526 +21 

1.6732 +20 

9.3870+16 

71.0000 

71.8019 

1.0425 +21 

- 

2.7966 +20 

1.2469 +19 

6.9957+15 

84.8520 

86,0000 

1.1298 +20 

8.6+16 

3.0309+19 

1.35140+18 

7.5817+14 
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Table 12 

Sound Speed, Dynamic Viscosity, Kinematic Viscosity, and Thermal Conductivity 
versus Geopotential Altitude, 0 to 86 km 


Altitude 

Sound 

Speed 

Dynamic Viscosity 

Kinematic Viscosity 

Thermal G 

nductivity 

H 

Z 

c s 

P 

0 

V 

9/9 o 

k 

k/k 0 

km' 

km 

m/s 

kg/(m • s) 


m 2 /s 


W/(m • K) 

- 

0.0000 

0.0000 

340.30 

1.7894 -5 

1 .0000 +0 

1.4607 -5 

l .0000 +0 

2.5326 -2 

1.0000 +0 

11.0000 

1 1 .0190 

295.07 

1.4216-5 

7.9447 - 1 

3.9064 5 

2.6743 +0 

1.9505-2 

7.7015-1 

20.0000 

20.Q631 

295.07 

1.4216-5 

7.9447 -1 

1.6148 -4 

I .1055 +1 

1 9505-2 

7.7015-1 

32.0000 

32.1619 

303.13 

1 .4868 - 5 

8.3090 - 1 

1 . 1 242 - 3 

7.6964+1 

2.0510-2 

8.0983 - 1 

47.0000 

47.3500 

329.80 

1.7037 -5 

9.5211-1 

1.1935 -2 

8.1703 +2 

2.3938 -2 

9.4521 -1 

51.0000 

51.4124 

329.80 

1.7037 -5 

9.5211 - 1 

1.9773-2 

1.3537 +3 

2.3938 -2 

9.4521 -1 

71,0000 

71 8019 

293.71 

1.4106 5 

7.8832-1 

2.1968-1 

1.5039 +4 

i. 9336 -2 

7.6349 - 1 

84.8520 

86.0000 

274.04 

1.2529-5 

7.0018-1 

1 8007 +0 

1 ,2327 +5 

1 .6962 - 2 

6.6976 - 1 


Table 1 3 

Temperature, Pressure, and Density versus Geometric Altitude, 86 to 1000 km 


Altitude 

Temperature 

Pressure 

Densi 

y 

Z 

H 

T 

T 

T « 

P 

P 

p / p o 

P 

p/p 0 

km 

km' 

K 

B c 

K 

mbar 

lorr 


kg/ m J 


86.0 

84.852 

186.87 

-86.28 

1 86.95 

3.7338- 3 

2.8006- 3 

3.6850 - 6 

6.958- 6 

5.680- 6 

91.0 

89.716 

186.87 

-86.28 

187.36 

1.5380- 3 

1.1536- 3 

1.5179- 6 

2.860- 6 

2.335 - 6 

110.0 

108.129 

240.00 

-33.15 

254.93 

7.1040- 5 

5.3284 - 5 

7.011 1 - 8 

9.708- 8 

7.925- 8 

120.0 

117.777 

360.00 

86.85 

397.91 

2.5382 - 5 

1.9038- 5 

2.5050- 8 

2.222- 8 

1.814- 8 

150.0 

146.542 

634.39 

361.24 

762.34 

4.5421 - 6 

3 .4069- 6 

4.4827 - 9 

2.076- 9 

1.694- 9 

200.0 

193.899 

854.56 

581.41 

1161.84 

8.4732- 7 

6.3554- 7 

8.3624-10 

2.541 -10 

2.074-10 

300.0 

286.480 

976.01 

702.86 

1594.82 

8.7699- 8 

6.5780- 8 

8.6552-11 

1.916-11 

1.564 -11 

400.0 

376.320 

995.83 

722.68 

1804.53 

1.4517- 8 

1.0888- 8 

1 .4327-11 

2.802-12 

2.288 - 1 2 

500.0 

463.539 

999.24 

726.09 

2019.70 

3.0234- 9 

2.2677 - 9 

2.9834-12 

5.215-13 

4.257-13 

600.0 

548.252 

999.85 

726.70 

2517.13 

8.2126-10 

6.1599-10 

8.1052-13 

1 137-13 

9.278-14 

700.0 

630.563 

999.97 

726.82 

3621.33 

3.1907-10 

2.3932-10 

3 1489- 13 

3.069-14 

2.506-14 

800.0 

710.574 

999.99 

726.84 

5225.13 

1.7036-10 

1.2778-10 

1.6813-13 

1.136-14 

9.272 -55 

900.0 

788.380 

1000.00 

726.85 

6577.15 

1.0873 -10 

8.1555 -11 

1.0731 -13 

5.759-15 

4.701-15 

1000.0 

864.071 

1000.00 

726.85 

7351.17 

7.5137-11 

5.6357-11 

7.4154-14 

3.561-15 

2.907 -16 
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Table 14 

Seven Secondary Atmospheric Properties versus Geometric Altitude, 86 to 1 000 km 


Altitude 

Acceleration 
Due to 
Gravity 

Pressure 

Scale 

Height 

Number 

Density 

Mean 

Particle 

Speed 

Mean 

Collision 

Frequency 

Mean Free 
Path 

Mean 

Molecular 

Weight 





Z 

H 

£ 

H p 

N 

V 

V 

L 

M 

km 

km' 

m/s 2 

kin 

m' 3 

m/s 

s' 1 

m 

kg/kmol 

86.0 

84.852 

9.5466 

5.621 +0 

1 .447 +20 

369.7 

3.17+4 

1.17 -2 

28.95 

91.0 

89.716 

9.5318 

5.642 +0 

5.962 +19 

370.1 

1.31+4 

2.83-2 

28.89 

1 10.0 

108.129 

9.4759 

7.723 +0 

2.144 +18 

431.7 

5.48 +2 

7.88 - 1 

27.27 

120.0 

117.777 

9.4466 

1.209+1 

5.107 +17 

539.3 

1.63+2 

3.31 +0 

26.20 

150 0 

146.542 

9.3597 

2.338 +1 

5.186+16 

746.5 

2.3 +1 

3.3 +1 

24.10 

200.0 

193.899 

9.2175 

3.618 +1 

7.182 +15 

921.6 

3.9 +0 

2.4 +2 

21.30 

300.0 

286.480 

8.9427 

5.119+1 

6.508 +14 

1079.7 

4.2 -1 

2.6 +3 

17.73 

400.0 

376.320 

8.6799 

5.968 +1 

1.056 +14 

1148.5 

7.2 -2 

1.6 +4 

15.93 

$00.0 

463.539 

8.4286 

6.879 +1 

2.192 +13 

1215.1 

1.6 -2 

7.7 +4 

14.33 

600.0 

548.252 

8.1880 

8.824 +1 

5.949 +12 

1356.4 

4.8 -3 

2.8 +5 

ll.Sl 

700.0 

630,563 

7.9576 

1.306+2 

! 2.311 +12 

1627.0 

2.2 -3 

7.3 +5 

8.00 

800.0 

710.574 

7.7368 

1.939+2 

1.234+12 

1954.3 

1.4 -3 

1 ,4 +6 

5.54 

900.0 

788.380 

7.5250 

2.509 +2 

7.876 + 11 

2192.7 

1.0 -3 

2.1 +6 

4.40 

1000.0 

864.071 

7.3218 

2.882 +2 

5 442 +11 

2318 1 

7.5 -4 

3.1 +6 

3.94 


Table 1 5 

Number Densities of Six Atmospheric Species versus Geometric Altitude, 86 to 1000 km 


Altitude 



Number Density 



Z 

H 

n(N 2 ) 

n(O) 

n(0 2 ) 

n(Ai) 

n (He) 

nOD 

km 

km’ 

in 3 

m' 3 

m -3 

m' 3 

m’ 3 

m 3 

86.0 

84.852 

1.130+20 

8.600+16 

3.031 +19 

1.351 +18 

7.582 +14 


91.0 

89,716 

4.643 +19 

2.953 +17 

1.234 +19 

5.478 +17 

3.419+14 


1 10.0 

108.129 

1.641 +18 

2.302 +17 

2.621 +17 

1 ,046 + 1 6 

5.821 +13 


120.0 

1 17,777 

3.726+17 

9.274+16 

4.395 +16 

1.366 +15 

3.888 +13 


150.0 

146.542 

3.124+16 

1.780+16 

2.750+15 

5.000 + 1 3 

2.106 +13 

3.767 +1 1 

200.0 

193.899 

2.925 +15 

4.050 +15 

1.918 +14 

1 .938 + 1 2 

1.310+13 

1.630+1 1 

300.0 

286.480 

9.593 +13 

5.433 +14 

3.942 +12 

1.568 +10 

7.566 +12 

1.049 +1 1 

400.0 

376.320 

4,669+12 

9.584 +13 

1.252 +11 

2.124+ 8 

4.868 +12 

8.961 +10 

500.0 

463.539 

2.592 +11 

1.836 +13 

4.607 + 9 

3 446 + 6 

3.215 +12 

8.000+10 

600.0 

548,252 

1.575 +10 

3.707 +12 

1.880+ 8 

6.351 + 4 

2.154 +12 

7.23! +10 

700.0 

630.563 

1.038+ 9 

7.840 +11 

8.410 + 6 

1.313+ 3 

1 .46 1 + 1 2 

6.556+10 

800.0 

710.574 

7.377 + 7 

1.732 +1 1 

4.105+ 5 

3.027 + 1 

1.001 +12 

5.961 +10 

900.0 

788.380 

5.641 + 6 

3.989 +10 

2,177 + 4 

7.742 - 1 

6.933 +1 1 

5.434 +10 

1000.0 

864,071 

4.626+ 5 

9.562 + 9 

1.251 + 3 

2.188- 2 

4.850+1 1 

4.067 +10 
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APPENDIX A 


BOUNDARY-VALUE NUMBER DENSITIES OF 
ATMOSPHERIC CONSTITUENTS 


The boundary-value neutral number densities of the several constituents defined to com- 
prise the U.S. Standard Atmosphere at 86 km and above were determined using a deductive 
process based upon several assumptions. It was decided to include as constituents of this 
model atmosphere only those species which are known to contribute significantly to the 
total number density in any portion of the atmosphere between 86 and 1000 km, because 
of either their mixing distribution below the turbopause or their diffusive distribution above 
this height. Those gases which appear never to contribute more than about 0.5 percent of 
the total composition at any height within this region, or which for various reasons do not 
exhibit predictable behavior, were purposely omitted. Using these guidelines, the following 
gases were included: molecular nitrogen, N 2 ; molecular oxygen, 0 2 ; argon, Ar; helium, He; 
and atomic oxygen, O. Atomic hydrogen, H, was included at heights above 149 km, but 
was not included in boundary-value considerations at 86 km. The remaining neutral gases 
which were used in establishing the sea-level value of the mean molecular weight, but which 
are not used at other heights in this model, are listed with the major gases and their respec- 
tive contribution to the sea-level mean molecular weight in table 16. 


Table 16 

Sea-level Atmospheric Composition 


Species 

Fractional 

Volume 

F. 

Molecular Weight 
of Species 

Mj 

F.M. 
1 1 

N 2 

0.78084 

28.0134 

21.87398326 

0 2 

0.209476 

31.9988 

6.70298063 

Ar 

0.00934 

39.948 

0.37311432 

CO, 

0.000314 

44.00995 

0.01381912 

Ne 

0.00001818 

20.183 

0.00036693 

He 

0.00000524 

4.0026 

0.00002097 

Kr 

0.00000114 

83.80 

0.00009553 

Xe 

0.000000087 

131.30 

0.00001142 

ch 4 

0.000002 

16.04303 

0.00003208 

H, 

0.0000005 

2.01594 

0.00000101 


2F. = 0.99999714 
1 


2(F. • = 28.964425 
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The first three of the gases used in this model comprise more than 0.9996 of the air in any 
unit volume at sea level, as is evident from summing the fractional composition F. over 
these three species in table 16. Since the fractional volumes of these major species do not 
change significantly below the mesopause, which in this model is located at 86-km altitude, 
the sea-level fractional composition can be assumed to be approximately correct at 86 km. 
It is believed, however, that photochemical processes lead to small quantities of atomic 
oxygen in this height region, and a fractional amount of about 0.00059 by volume, or 
exactly 8.6 X 10 16 atoms per m 3 , was agreed upon by Task Group II as an acceptable con- 
centration of O for 86 km. 

The introduction of the fixed amount of atomic oxygen at this height, and the simultane- 
ous elimination of some minor species made it necessary to adjust the fractional concen- 
trations of each of the four remaining species from their known sea-level values, F., by a 
common unknown factor, e, to the 86-km fractional composition values, F.' , such that 

F! = e * F { (58) 

The 86-km fractional composition of atomic oxygen is equal to the ratio of n(O) to N, 
where n(O) is equal to 8.6 X 10 16 m' 3 , the adopted atomic-oxygen number density for 
that height, and N is the unknown total number density at that height. The sum of the 
86-km fractional composition of the remaining four species is 

F i (59) 

4 4 

The sum of the total of the 86-km fractional compositions, (that is, of the Five species 
adopted to comprise the model at this height) must equal unity in accordance with the 
expression: 


6 



n(Q) 

N 


(60) 


The total number density, N, is expressible in terms of the mean molecular weight, M, 
Avogadro’s constant, N. , and the mass density, p, the value of which is known at 86 km 
from other considerations. This relationship is 


N = 



(61) 


The mean molecular weight at 86 km is the sum of the products F| • M. over the five gases 
comprising the model at this height. For the atomic oxygen this product is 


F'(0)-M(0) 


n(0)-M(0) 

N 


(62) 
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while the sum of the products of the remaining four gases is expressible as 



such that the mean molecular weight, M, at 86 km is expressed as 


M = 



F.-M; 


n(0>M(0) 

_ 


( 63 ) 


(64) 


Eliminating M between equations 61 and 64, and solving for e yields 


e 


N a -p - n (0) • M (0) 



(65) 


The elimination of e between equations 60 and 65 leads to the following expression for 
total number density: 


N 



+ n(0) 


( 66 ) 


From table 16 the value of 2F., the sum of F. for the four species N 2 , 0 2 , Ar, and He is 
seen to be 0.99966124, while the value of 2F. M. for the same four species is seen to be 

4 l I 

28.95009918. The value of M(O) is taken to be one half of the value of M(0 2 ), also given 
in table 16. The value of 8.6 X 10 16 m' 3 was adopted for n(O), as previously stated, and 
N a has the standard value 6.022169 X 10 26 kmol -1 . The value of p at 86 km is found to 
be 6.957880 X 10‘ 6 kg/m 3 . These values, introduced into equation 66, yield a number 
density of 1.447265 X 10 20 m' 3 at 86 km. This value introduced into equation 58 leads 
to e = 0.99974445, while equation 64 then yields M = 28.952208 for the molecular weight 
at 86 km. 


For the 86-km height, the values of F' the fractional composition of each of the five spe- 
cies comprising the model at that height, are given in table 1 7 as the product e • F[, along 
with the corresponding products F| • M p and the corresponding values of N • F', the num- 
ber densities of the five gas species comprising the model atmosphere at 86-km height. 

The value of 2F.', the sum of the Five values of F.' listed in table 17 is seen to be 0.999999999, 
essentially the unity value which it should have. The sum of the five values of F. • \F and 
of the five values of n ; , that is, SF.' • M.' and 2n., both of which are also given in table 17, 
shows essentially exact agreement with the value of their respective equivalents, M and N, 
computed independently. Thus, the validity of the computation is established. 
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Table 1 7 

Number Densities and Molecular Weight at 86 km 


■ 

F! = eF. 
1 1 

F.' * m. (kg/kmol) 

n. = F[ • N (m' 3 ) 

N 2 

0.7806404557 

21.86839334 

1.129793736 X 10 20 

0 2 

0.2094224682 

6.701267675 

0.303 0898426 X 10 20 

Ar 

0.00933761315 

0.3730189704 

0.0135140022 X 10 20 

He 

0.00000523866 

0.0000209683 

0.0000075817 X 10 20 

O 

0.00059422421 

0.0095072308 

0.00086 X 10 20 


2 F.' = 0.99999999992 

2Ff • M. = 28.9522082 

2n. = 1.447265163 X 10 20 


5 ‘ 

5 1 ‘ 

5 1 












APPENDIX B 


A SEGMENT OF AN ELLIPSE TO EXPRESS 
TEMPERATURE VERSUS HEIGHT 


It is desired to determine the expression for a temperature function for a limited height 
region, Z g = 91 to Z 9 = 1 10 km, in the plane defined by Z and T, such that the slope of 
the function at each of the end points exactly matches a prescribed value. At Z = Z g , 
where T = T g = 186.8673 K, the derivative of the function with respect to Z must be zero, 
to match the slope of the temperature-height profile in the isothermal layer between 86 
and 91 km. At Z = Z g , where T = T g = 240 K, the derivative of T with respect to Z must 
be 12 K/km to match the slope of a layer of constant temperature-height gradient between 
1 10 and 120 km. A suitably adjusted ellipse will satisfy these conditions. 

The general equation of an ellipse in terms of Z and T with center at Z = 0 and T = 0 is 


2 ?_ j2 

a 2 + A 2 


With the center shifted to Z = Z c and T = T c the expression becomes 

(Z - Z ) 2 (T - T ) 2 
1- + L_ = 1 

a 2 A 2 


(67) 


( 68 ) 


The derivative of equation 68 with respect to Z is 


2(2- Z e ) 2(T-T c ) dT 
a 2 A 2 


(69) 


To meet the condition for dT/dZ = 0 at Z = Z g , equation 69 is evaluated for those condi- 
tions, and it is found that Z c = Z g , such that equation 68 may be rewritten as 


(Z “ Z 8 ) 2 _ (T - T, ) 2 

a 2 

Evaluating equation 70 for Z = Z. and T = T. 

o o 

a = t 8 

ft 

Substituting Z c for its equal Z g in equation 69, and evaluating that expression for Z = Z g , 
where T = T g and where (dT/dZ) has the particular value L R g = 12, and finally solving 
the resulting expression for 1/a 2 yields 


A2 

leads to 

-T (71) 
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( 72 ) 



Since it was shown that Z = Z fi , the ellipse which meets the required derivative and temp- 

C o 

erature conditions, has its center at Z = 91 km and T = 263.1905 K, and equation 76 repre- 
sents the function which meets the required conditions. 
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APPENDIX C 


GEOPOTENTIAL AND MOLECULAR-SCALE TEMPERATURE 


The concept of geopotential as a measure of height was introduced by Bjerknes et al. (1910) 
who, according to Harrison (1954), made use of the term dynamic height in referring to 
the geopotential of a point because in meteorology, geopotential is preferable to geometric 
height as a representation of the vertical coordinate of the point. Bjerknes proposed the 
term geodynamic meter (gdm), or dynamic meter for short, as the name of the unit of 
dynamic height. With a slight revision in the definition of this concept, the Aerological 
Commission of the International Meteorological Organization in 1947 adopted the geo- 
potential meter (gpm) defined such that 1 gpm = 0.98 gdm = 9.8 m 2 /s 2 . This definition 
related the geopotential meter directly to the sea-level value of the acceleration of gravity. 

In the U.S. Standard Atmosphere of 1955, published in NACA Report 1235 (1955), the 
standard geopotential meter (sgpm) was defined in terms of the so-called standard sea-level 
value of the acceleration of gravity, 9.80665 m/s 2 , such that 1 sgpm = 9.80665 m 2 /s 2 . 

The use of the standard geopotential meter (now designated as m') was continued in the 
Air Research and Development Command (ARDC) Model Atmosphere, 1956 (Minzner 
and Ripley, 1956); in the U.S. Extension to the ICAO Standard Atmosphere (Minzner et al., 
1958); in the ARDC Model Atmosphere, 1959 (Minzner et al., 1959), in the U.S. Standard 
Atmosphere, 1962 (COESA, 1962); as well as in the 1975 revision of that Standard Atmo- 
sphere (COESA, 1975). Implicitly, geopotential has been used as a measure of height in 
all standard and model atmospheres, in which the tabulated values of atmospheric proper- 
ties have been calculated on the basis of a value of the acceleration of gravity which is in- 
variant with height. 

The numerical value of the height of a particular point above sea level, in regions where 
the sea-level value of the acceleration of gravity is equal to or less than 9.80665 m/s 2 , is 
smaller when expressed in geopotential than when expressed in geometric meters. One 
standard geopotential meter is exactly equal to one geometric meter only under very spec- 
ialized and perhaps hypothetical conditions. These conditions involve the existence of a 
region with a gravitational-field value of exactly 9.80665 m/s 2 (generally near the earth’s 
surface), and with a zero vertical gradient of that field over a height interval of at least one 
geometric meter. The concept of 1 m' is more precisely defined, however, in terms of a 
specific change in potential energy of a unit mass, such that 1 m' is exactly the height 
increment through which one must lift one kilogram (mass) in order to increase its poten- 
tial energy by 9.80665 J. The geometric equivalent of this height increment is essentially 
one geometric meter (m) at sea level, but increases to more than 1 m with increasing height, 
in a manner inversely proportional to the corresponding decrement of the acceleration of 
gravity. 
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The concept of geopotential has some very meaningful uses in meteorology and aeronomy. 
One of its principal advantages is the reduction from 4 to 3 of the number of variables in 
the hydrostatic equation when geopotential is introduced. This feature is particularly de- 
sirable in expressions involving the integrated forms of the hydrostatic equation. Thus, the 
use of this concept, either implicitly or explicitly in the definition of standard atmospheres 
produced prior to the development of high-speed digital computers, simplified the defining 
equations and reduced the work involved in manual computation of the related tables. 

The use of geopotential has been maintained in all subsequent standard atmospheres in 
order to avoid any revision of the lower portion of the tables which have represented the 
established standard for the past 50 years. 

The introduction of the concept of molecular-scale temperature, T M , came at a time when 
standard-atmosphere tables were being extended to heights where the composition, and 
hence the mean molecular weight, were unknown. The concept (without a particular name) 
was first applied to atmospheric models by the Rocket Panel (1952) in a paper which re- 
ferred to a “derived ‘Temperature’ based upon an assumed constant mean molecular weight 
H of 28.966 g/mole for the atmosphere.” The use of this special kind of temperature not 
only avoided the problem of determining or guessing at a value of mean molecular weight 
at high altitudes, but also eased the problem of hand or desk-computer calculations by 
leading to simpler equations than would have resulted if specific height functions had been 
introduced for both kinetic temperature and mean molecular weight. 

The name molecular-scale temperature for this derived temperature first appeared in the 
ARDC Model Atmosphere, 1959. In this model, molecular-scale temperature, T M , at any 
height Z was related to the kinetic temperature, T, at height Z by the relationship T M (Z) 

= T(Z) [M 0 /M(Z)] where M 0 was the currently accepted sea-level value of the mean mole- 
cular weight of air, 28.9644 kg/kmol, and M was the implicit value of the mean molecular 
weight of air at height Z. This concept was carried into the U.S. Extension to the ICAO 
Standard Atmosphere, 1958; the ARDC Model Atmosphere, 1959; the U.S. Standard Atmo- 
sphere, 1962 (COESA, 1962), and has been carried into the 1975 revision to that Standard 
Atmosphere (COESA, 1975). 

The temperature-height profile of each of the first three of these four previously published 
models was defined in terms of T M (in units of Kelvin, K) as a function of geopotential 
height H (in units of m') over the entire height range of these models. Since many aerono- 
mers studying the upper atmosphere are unfamiliar with these quantities and prefer the 
quantities kinetic temperature, T, and geometric heights, Z, it was decided that, in the 
U.S. Standard Atmosphere, 1962 (COESA, 1962) the model should be divided at 90 geo- 
metric kilometers, with temperatures in that part of the model from 0 to 90 km defined in 
terms of T M and H, and in that part above 90 km defined in terms of T and Z. This proce- 
dure has been continued in the 1975 U.S. Standard Atmosphere, but in this revision, the 
division occurs at 86 geometric kilometers where the geopotential is 84.8520 km', and 
M = 28.9522 kg/kmol. 
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Before the name molecular-scale temperature was used in model atmospheres, Brombacher 
(1953), in an attempt to simplify the calculation of standard atmospheres, combined the 
height-dependent acceleration of gravity, the height-dependent molecular weight, and the 
height-dependent kinetic temperature into a single variable, which he called scale-height 
temperature. The use of this variable, however meritorious, was never adopted in U.S. 
Standard Atmospheres. 
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APPENDIX D 


THE CALCULATION OF A DYNAMIC MODEL FOR THE 
1975 U.S. STANDARD ATMOSPHERE 


INTRODUCTION 

The objective of this appendix is to describe the procedure for the calculation of a dynamic 
model of the earth’s atmosphere between 50 and 1 50 km, comprised of an internally con- 
sistent set of diumally averaged gas concentrations versus altitude, one for each of the four 
major atmospheric gas species in this height region, that is, nitrogen, molecular and atomic 
oxygen, and argon, where each of these concentration profiles meets the following two 
conditions: 

1. The concentration values versus height are the result of a time-dependent, photo- 
chemical-transport calculation which incorporates measured chemical-reaction-rate 
constants, solar-radiation fluxes, and turbulent diffusion coefficients into coupled 
sets of equations of motion and continuity. 

2. The calculated number densities of each of the four species at 1 50 km fall within 
particular limits recommended by the COESA Working Group. 

This sophisticated and detailed calculation serves to establish the physical basis for the 
generation of dynamic models of the earth’s atmosphere, and yields height profiles of 
number-density flux values which are approximated by artificially adjusted functions for 
the calculation of the 1975 U.S. Standard Atmosphere. 

BASIC CONSIDERATIONS 

The species considered are O, 0 2 , 0 3 , O 1 D, 0 2 ( l A g ), OH, H, H0 2 , H 2 0, H 2 0 2 , H 2 , Ar, 
and He. The number densities, from 50 to 1 50 km, are obtained through a semi-implicit, 
finite-difference solution of a system of mass- and momentum-conservation equations 
(Shimazaki, 1967; Keneshea and Zimmerman, 1970). In these calculations, thermal-diffu- 
sion factors for the species, H, H 2 and He have the values respectively of -0.39, -0.31, and 
-0.36 (Zimmerman and Keneshea, 1975). The numerical approach is essentially that in- 
troduced by Shimazaki (1967), but modified at the boundaries and in the volume integra- 
tions, following George et al. (1972). Table 18 lists the chemical reactions and the asso- 
ciated rate constants actually used in the generation of the resulting concentration profiles. 
It should be noted, however, that refined rate-constant measurements, made since these 
calculations were completed, indicate the need for a revision of some of the listed values. 
The current calculations have not been updated with these new rate-constant values, how- 
ever, since the changes have only a negligible influence on the concentrations of O, 0 2 , Ar, 
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Table 18 

The Chemical Reactions and Associated Reaction Rates kj 
Expressed in the Form of the Value of Rate Coefficients 
A j5 B., and Cj where k. = A. • (T/300) B i X exp(Cj/T) 


REACTION 

A i 

B , 

c. 

J 

1 

O 

4 O 

4 M 

-o 2 

4 M 

3.00E-33 

-2.9 


2 

O 

+ o 2 

4 M 

-°3 

4 M 

5. 50 E- 34 

-2.6 


3 

o 

+ °3 


•°2 

+ 0 2 

1.20E-11 


-2.00E403 

4 

H 

+ o 3 


-o 2 

4 OH 

2.60E-11 



5 

OH 

+ O 


- H 

+ °2 

5.00E-11 



6 

OH 

+ o 3 


- ho 2 

4 0 2 

4.00E-14 



7 

H 

+ o 2 

4 M 

- ho 2 

4 M 

7.40E-33 


6.10E+02 

8 

HO 2 

4 O 


- OH 

4 O, 

1.00E-11 



9 

h° 2 

+ ° 3 


- OH 

+ °2 + °2 

1.00E-17 



10 

OH 

4 OH 


- h 2 o 

4 O 

2.00E-12 



11 

OH 

+ ho 2 


* h 2 o 

+ o 2 

2.00E-10 



12 

H 

+ ho 2 


- H 2 

+ 0 2 

3.00E-12 



13 

H 

+ ho 2 


- OH 

4 OH 

1.00E-11 



14 

O 

+ H 2 


- OH 

4 H 

7.00E-11 


-5.10E+03 

15 

ho 2 

+ ho 2 


- H 2°2 

* °2 

3.00E-12 



16 

OH 

+ h 2 o 2 


- h 2 o 

4 H0 2 

1.70E-11 


-9.00E+02 

17 

O 

+ h 2 o 2 


- OH 

4 H0 2 

4.00E-15 



18 

H 

+ h 2 o 2 


- H 2 

4 H0 2 

3.90E-11 


-4.60E+03 

19 

O 1 D 

+ ° 3 


- o 2 

4 0 2 

3.00E-10 



20 

O 1 D 

+ °2 


- o 

+ °2 

6.00E-11 



21 

O 1 D 

+ n 2 


- o 

+ N 2 

9.00E-11 



22 

0*D 

+ H 2 


- OH 

4 H 

1.00E-11 



23 

O 1 D 

+ HjO 


- OH 

4 OH 

1.00E-11 



24 

0 2 l Ag + 0 3 


-°2 

+ 0 3 

3.00E-15 



25 

0 2 j A g + M 


-°2 

4 M 

4.40E-19 



26 

Oj'Ag + H 


- OH 

4 O 

1.10E-14 



27 

o 2 1 a 

g 


- o 2 


2.58E-04 



28 

°2 

+ hv 


- o 

4 O 




29 

o 2 

+ hi' 


— O 1 D 

4 0 




30 

o 3 

4 hi; 


-° 2 

4 O 




31 

o 3 

+ hi; 


- O 1 D 

4 0 2 1 A g 




32 

h 2 o 

+ hv 


- OH 

4 H 




33 

h 2 o 2 

4 hi' 


- OH 

+ OH 





Note: The units of the two- body reaction rates are cm 3 s' 1 , while those 
for the three-body reaction rates are cm 6 s' 1 . 




and He in the altitude region above 80 km. The intensity of the solar-radiation flux used 
in these calculations is 0.65 of that shown in figure 17, which depicts the Ackerman (1971) 
values of solar-radiation flux versus wavelength. The absorption cross sections were taken 
from various sources. For 0 2 and 0 3 , these cross sections were taken from the compila- 
tion of Ackerman (1971) with the exception of those for the Schumann- Runge bands of 
0 2 , for which region the values measured by Hudson and Mahle (1972) were used. The 
adopted absorption cross sections for water vapor and hydrogen peroxide are those re- 
ported by Watanabe and Zelikoff (1953), and by Volman (1963), respectively. The temp- 
erature-height profile up to 1 50 km and the values of mean molecular weight up to the 
turbopause are those recommended by the Working Group of COESA. Using these data, 
the initial species distributions were calculated assuming complete mixing up to the turbo- 
pause, and diffusive equilibrium above it. 


SOLAR FLUX 



Figure 17. Solar radiation flux versus wavelength in the region 
from 1 1 5 to 31 0 nm as reported by Ackerman. 

The total number density was obtained by integrating the hydrostatic equation, where the 
sea-level values of mass density and of mean molecular weight were taken from the U.S, 
Standard Atmosphere , 1962 (COESA, 1962). 

Beginning with these static profiles, the steady-state solution of all species was determined. 
The time-dependent calculations were then allowed to proceed for 1 5 solution days using 
a semi-implicit, finite-difference technique, a variable-time step up to 30 minutes, and a 
Fixed-height step of 100 m. This stringent height step was shown to be necessary to restrict 
the errors generated by species gradients when height steps larger than 100 m were used. 
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The height-dependent, turbulent diffusion coefficients used are shown in figure 18, and 
are based upon observations of turbulence in chemical trails (Philbrick et al., 1973). These 
values are derived (Zimmerman and Trowbridge, 1973) from the fluctuation dynamics ob- 
served in rocket-borne chemical releases, and are valid from about 88 to 112 km. Because 
of the lack of chemical-tracer wind and turbulence measurements in the altitude region 
between 50 and 88 km, an exponential fit has been assumed between the reported value 
of 1 X 10 5 cm 2 s' 1 (Beaudoin et al., 1967) at 50 km, and the values at 88 km. 



EDDY - DIFFUSION COEFFICIENT, m 2 /s 

Figure 18. Eddy-diffusion coefficient versus altitude. 

RESULTS 

The time-dependent calculations were continued for the above mentioned period of time, 
after which the species concentrations reproduced themselves to within 1 percent over a 
diurnal cycle, a condition which is called arriving at diurnal reproducibility. The diurnal 
averages of the concentration of O, 0 2 , and Ar are then calculated and extrapolated to 
250 km by assuming diffusive equilibrium without thermal diffusion above the 1 50-km 
boundary. Figure 19, depicting the height profiles of the N 2 concentration and tempera- 
ture, shows the initial conditions used in these one-dimensional calculations. Figure 20 
shows the resulting diumally averaged height profiles of 0, 0 2 , and Ar, each of which is in 
good agreement with the 1 50-km values recommended by the COESA Working Group, and 
shown as error bars. 
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N 2 CONCENTRATION, m ' 3 

10 14 10 15 10 16 10 17 10 18 10 19 10 2C 10 21 10 22 10 23 



Figure 19. Altitude profile of kinetic temperature and molecular nitrogen concentration 
with the range of 150-km values recommended by the COESA Working Group indicated 
by an error bar. 



Figure 20. Altitude profiles of diurnally-averaged concentrations of O, 0 2 , and Ar, 
with the range of 150-km values recommended by the COESA Working Group indicated 
by corresponding error bars. 
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Thus, it has been demonstrated that an internally self consistent model of the density struc- 
ture of the upper mesosphere and lower thermosphere may be calculated from measured 

values of solar radiation flux, chemical-reaction-rate constants, and derivatives of measured 

vertical-turbulent-transport parameters deduced from chemical-trail studies. 
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